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Abstract: We demonstrate that there exists a large class of J-{R,G) action functionals 
of the scalar curvature and of the Gaufi-Bonnet invariant which are able to relax dynam- 
ically a large cosmological constant (CC), whatever it be its starting value in the early 
universe. Hence, it is possible to understand, without fine-tuning, the very small current 
value Aq ~ Hq of the CC as compared to its theoretically expected large value in quantum 
field theory and string theory. In our framework, this relaxation appears as a pure grav- 
itational effect, where no ad hoc scalar fields are needed. The action involves a positive 
power of a characteristic mass parameter, Ai, whose value can be, interestingly enough, 
of the order of a typical particle physics mass of the Standard Model of the strong and 
electroweak interactions or extensions thereof, including the neutrino mass. The model 
universe emerging from this scenario (the "Relaxed Universe") falls within the class of 
the so-called AXCDM models of the cosmic evolution. Therefore, there is a "cosmon" 
entity X (represented by an effective object, not a field), which in this case is gener- 
ated by the effective functional J-{R, Q) and is responsible for the dynamical adjustment 
of the cosmological constant. This model universe successfully mimics the essential past 
epochs of the standard (or "concordance") cosmological model (ACDM). Furthermore, 
it provides interesting clues to the coincidence problem and it may even connect naturally 
with primordial inflation. 
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1. Introduction 

The status of cosmology as a quantitative subdiscipline of physics is largely based on 
the high-quality observations made during the last two decades. The availability of these 
data enabled cosmologists to demonstrate that cosmological models are able to explain 
not only qualitative features of the universe and its dynamics, but also quantitatively 
reproduce the details of the many observed cosmic phenomena [1]. Yet, the discovery 
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of the accelerated expansion of the universe at the present epoch [2] revealed a gap in 
our understanding of the universe. The search for the underlying mechanism producing 
the accelerated expansion revived an old problem of theoretical physics: the cosmological 
constant problem [3, 4], which manifests itself as the double conundrum on the tiny observed 
value of the cosmological constant (CC) in Einstein's equations - the "old CC problem" [5] 
- and the puzzling fact that this value is so close to the current matter density - the 
"cosmic coincidence problem" [6]. Even though the simplest candidate for the acceleration 
mechanism is the existence of a small positive cosmological constant, A, now incorporated 
into the benchmark model of modern cosmology, or "concordance" ACDM model, there 
is quite a list of challenging cosmological problems of various sorts afflicting the standard 
model of cosmology, see e.g. [7]. 

Above all these problems there is a problem of greatest importance, one which is by 
far the most troublesome one from the point of view of Fundamental Physics, to wit: 
the size of the cosmological constant energy density required to describe the cosmological 
data is in drastic disaccord with the value predicted by any known fundamental physical 
theory based on quantum field theory (QFT) or string theory. Indeed, from a theoreti- 
cal point of view, the cosmological constant energy density, pA = A/(87rG7v) {Gn being 
Newton's constant), is a single number, but of remarkable composition as it involves all 
sources of vacuum energy: e.g. the zero-point energy density of a quantum field of mass 
m produces a contribution to the CC energy density of the order of m^, where bosonic 
fields contribute positively and fermionic fields negatively; phase transitions in the early 
universe also leave an imprint on the value of pA in the form of vacuum energy associated 
to spontaneous symmetry breaking (e.g. the energy density of the electroweak vacuum); 
or the non-perturbative condensates (as e.g. the quark and gluon condensates of the QCD 
vacuum) etc. All these contributions, which differ in size and sign, have a common char- 
acteristic: they are all much larger in absolute value than the experimentally determined 
vacuum energy density, ~ (2.3 x lO^^eV)^ [1]. There is a (remote) possibility that 
these manyfold effects might cancel among themselves to produce the observed value of 
the CC. This, however, requires extreme fine-tuning which, in perturbative QFT, needs 
to be iterated to all orders of the perturbative expansion. Such explanation, therefore, 
although technically possible, is utterly unconvincing. 

If we discard the explanations based on fine-tuning, and we take seriously the idea that 
all the contributions to the vacuum energy do have an impact on the final CC value, then 
we are unavoidably left with a value of p\ many orders of magnitude larger than p^. Since 
we have no experimental verification of the dynamics in the UV limit, the size and sign of 
the "initial" CC energy density - which we may call p\ - are not precisely known, except 
the fact that its absolute value must be really huge when measured in units of the observed 
CC energy density; i.e. the number |paI/Pa ™ust be enormous. 

To reconcile the anomalously large size of this ratio with the behavior of the observed 
universe, the value of p\ has to be neutralized in some reasonable manner. One possibility 
would be the existence of a mildly broken symmetry requiring the CC to be very small 
(zero if the symmetry would be exact). Such a symmetry, however, is not presently known. 
Supersymmetry (SUSY) [8] , for instance, despite the first natural hopes [9] , cannot cure 
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the CC problem in any obvious way, because we know that SUSY must be broken at the 
electroweak scale or higher, which means that the SUSY particles must have masses larger 
than their conventional counterparts, and as a result the "residual" vacuum energy left in 
the universe must again be of the order (or even larger) than that of the Standard Model 
(SM) of particle physics, i.e. of order M^, where M^r ~ 100 GeV is the scale of the weak 
gauge boson masses ^. Therefore, in general we expect, roughly, \pW/p\ ~ 10^^, which is an 
appallingly large number! Obviously, a theoretical breakthrough is demanded to account 
for this situation. 

From a physical point of view, it is especially appealing to have a dynamical adjustment 
of the CC value, because we can then avoid resorting to too preposterous a fine-tuning of 
the parameters ~ at least 55 digits in the SM case. This idea was originally pursued in 
terms of dynamical scalar fields [12], but it was later shown by S. Weinberg that it is 
generally obstructed by a "no-go theorem" [3]. Subsequently, the more modest idea of 
quintessence [13] was proposed without attempting to explain the smallness of the CC, but 
only to cope with some aspects of the cosmic coincidence problem. 

In this paper, we wish to face a possible dynamical solution that can escape the no- 
go theorem and that is based on a different concept, one that makes no use of scalar fields. 
Specifically, we refer to the concept of "dynamical effective A" or "effective vacuum energy" 
PAeff) in which, rather than replacing A by a collection of ad hoc ersatz fields, we stick to the 
idea that the CC term in Einstein's equations is still a "true cosmological term" , although 
we permit that "the observable CC at each epoch" can be an effective quantity evolving 
with the expansion of the universe: pAcff = PAcsiH), where H is the expansion rate or 
Hubble function. This general dynamical A approach has been recently re-emphasized in 
[14, 15], and it was explored with different particular implementations in the past from the 
point of view of QFT in curved-space time [16, 17], including some recent applications to 
structure formation [18] - see e.g. the review [19] and references therein. 

The variable CC approach was also tried long ago on purely phenomenological terms, 
cf. [20, 21]. Furthermore, the idea of a slowly running CC is in general compatible with the 
experimental data [22, 23], and has been even exploited as an alternative description of the 
notion of dynamical dark energy (DE) for a possible viable solution of the aforementioned 
cosmic coincidence problem [24, 25, 26]. But, most important of all, it may also be advan- 
tageously utilized as a powerful mechanism capable of solving (or highly alleviating) the 
big or "old" CC problem [3], viz. the absolutely formidable task of trying to explain the 
value itself of the current vacuum energy - not just its time evolution in the vicinity of 
it - on the face of its enormous input value p\ left over in the early times. 

Can we explain the measured value of A and at the same time insure that it remains 
safely small throughout the entire history of the universe? In this paper, we shall present 
a thorough attempt of this sort, actually one which substantially improves previous recent 

^In certain local realizations of SUSY, such as in Supergravity theories emerging from Eg x Es super- 
string theory with gluino condensation, one can have vanishing vacuum energy even if SUSY is broken [10]. 
However, this does not explain the small value of the CC, specially after one realizes of the vast "landscape" 
of possible string vacua [11], suggesting that there is no particular reason for a given vacuum choice in this 
framework. 
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attempts along these lines [27, 28]. From now on, we shall call this new process of dynamical 
neutralization of the "true" cosmological constant: the "relaxation mechanism" of the CC. 

Relaxation offers indeed an entirely different perspective on the old CC problem. Many 
CC adjustment mechanisms, partly discussed above, function unilaterally, as e.g. when 
using scalar fields. In other words, in these cases the dynamics of the relaxation mechanism 
works "intrinsically" to neutralize the real value of the CC. A new quality could be added 
to the relaxation mechanism if it becomes bilateral. Namely, apart from the mechanism 
itself for neutralizing the value of the CC, we could also envisage that the expansion of 
the universe, fueled by the big CC value p^, enables the very action of the relaxation 
mechanism. The idea is that the relaxation mechanism neutralizes the value of but at 
the same time the existence of this value sets the relaxation mechanism to motion. Since 
the expansion of the universe should be responsible for the triggering of the relaxation 
dynamics, it is reasonable to assume that such mechanism should be closely related to the 
expansion itself, or more precisely, to the interaction governing the expansion, i.e. gravity. 

Despite the appeal of such a feedback, a number of important questions have to be 
answered: i) Is it possible that the observed value of the CC (or effective vacuum energy) is 
small exactly because the "real" value of the CC is large? ii) Is there a way to dynamically 
counterbalance the effects of a large CC in all epochs of the cosmic evolution? iii) Can the 
feedback between the relaxation mechanism and expansion be formulated in terms of an 
action functional? In this paper, we present a modified gravity theory based on an action 
functional J-{R,G) of the Ricci scalar, R, and the Gauss-Bonnet invariant, Q, in which all 
the above questions have an affirmative answer. 

A simple model exploring the viability of a mechanism cancelling dynamically the ef- 
fects of a large CC was presented in [27]. Previously, it had been shown that the inho- 
mogeneous equation of state (EOS) of a cosmic fluid can be interpreted as an effective 
description of modified gravity theories or even braneworld models [29]. In a cosmological 
model containing only two components, a large cosmological constant of an arbitrary sign 
and a component with an inhomogeneous EOS of the form p = wp — I3H~°', for a > the 
cosmic expansion asymptotically tends to de Sitter phase with a small positive effective 
CC [27]. Essentially, the effect of a very large A is counterbalanced by the terms of the 
form when H is sufficiently small. In this cosmological model, at late times the 

universe expands as if there existed a small positive CC even though the actual CC is 
very large and it can even be negative. As there is no fine-tuning in this model, such a 
dynamical outcome can be rightfully identified as a viable approach for solving the cosmo- 
logical constant problem. One specific feature of the model deserves a special emphasis. 
Within this relaxation mechanism, the size of the CC is not only a problem, but also a 
part of the solution. To illustrate this fact let us consider the asymptotic value of the 
Hubble parameter .ffasym ~ V^cff for a = 2. For this value of a we have -ffasym ~ 1/IPaI- 
The bigger p\ is, the smaller is the effective observed Acff. For some other recent work 
on relaxation mechanisms, see also [30] and [31]. For recent alternative ideas on the CC 
problem, see [32, 33]. 

The mechanism based on the inhomogeneous EOS, despite exhibiting CC relaxation 
properties, is incomplete since it does not contain matter or reproduce the observed se- 



-4- 



quence of epochs in the history of the universe. The idea of the CC relaxation mechanism 
has to be embedded into a reahstic cosmological model. Such a realistic implementation was 
made in [28] in the formalism of variable A interacting with a dark matter component, i.e. 
within the context of the so-called AXCDM model [24]. The variable CC contains contri- 
butions in the form of functions of curvature invariants such as R"^, W^^Rab and R"''"^'^ R^hcd- 
The model successfully produces the sequence of radiation, matter and de Sitter phases 
in the expansion of the universe. During the de Sitter phase, the action of a large CC is 
counterbalanced by the term l/H"^^ when H is sufficiently small, similarly as with having 
inhomogeneous EOS. However, during the matter and radiation dominated phases, the 
effect of a large A is equilibrated by the terms ~ — 1/2) and ~ l/(q — 1), respectively, 
where q is the deceleration parameter. Within this model one can easily produce an expan- 
sion history of the universe close to that of the ACDM model despite the persistence of the 
huge primeval CC, which is permanently subdued by the dynamical mechanism. Besides, 
a dedicated analysis of the growth of perturbations shows that the formation of structures 
in this kind of models is also comparable to the ACDM [34]. 

The model introduced and elaborated in [28] is thus phenomenologically very success- 
ful, but it lacks the action principle formulation. To fill this remaining gap, we proposed 
recently a specific model universe [35] based on an action functional for modified gravity. 
In the present paper, we extend and generalize this approach while maintaining the phe- 
nomenologically pleasing features discussed above. Specifically, we discuss here a modified 
gravity theory ("The Relaxed Universe") based on a class of J-{R,Q) action functionals 
capable of relaxing the highly "stressed" primeval state of our universe, namely a state 
which is beset by a large vacuum energy which prevents the startup of the normal ther- 
mal history. The relaxation mechanism allows not only to unblock this situation, it also 
operates very actively during all the subsequent epochs of the cosmic expansion without 
fine-tuning the parameters of the theory at any time. Furthermore, it involves a character- 
istic mass scale, Ai, whose value is of the order of a typical particle physics scale whether of 
the SM or of a Grand Unified Theory (GUT), and therefore avoids introducing extremely 
tiny masses (m^ Hq W~^^ eV) plaguing most proposals, as e.g. quintessence [13]. 

It should finally be emphasized that the cosmological model introduced here goes be- 
yond the idea of modified gravity with late time cosmic acceleration [36, 37]. Indeed, let 
us stress once more that the proposed gravity modifications inherent to our mechanism 
are crucial during the entire cosmic history, and not just in the late universe. The rea- 
son is that they primarily counterbalance the effect of a huge CC during all accessible 
epochs (radiation and matter), and without any fine-tuning. The late time acceleration 
phenomenon appears here only as a special, but certainly very important effect, either in 
quintessence-like, de Sitter or phantom-like disguise. 

The paper is organized as follows. In the next section, we sketch the old CC problem 
and its relation with the notion of fine-tuning. In section 3 we define and solve the J-{R, Q)- 
cosmology. In section 4 we discuss the working principle of the relaxation mechanism. 
The analysis of some concrete implementations of the model is performed in section 5, 
including detailed numerical results. In the last section we deliver our conclusions. Finally, 
in appendix A we provide information on our notation and conventions, in appendix B 
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we extend the discussion of the severity of the CC fine-tuning problem in QFT, and in 
appendix C we briefly discuss the relaxation mechanism in the solar system. 



2. The old CC problem as a fine-tuning problem 

Before we present our relaxation mechanism, let us summarize the old CC problem and 
discuss why the fine-tuning problem unavoidably appears if one does not modify some as- 
pects of the theory of gravity in interaction with matter. We wish to illustrate the problem 
within the context of the standard model (SM) of particle physics, and more specifically 
within the Glashow-Weinberg-Salam model of electroweak interactions. This is the most 
successful QFT we have at present (together with the QCD theory of strong interactions), 
both theoretically and phenomenologically, and therefore it is the ideal scenario where to 
formulate the origin of the problem. As is well-known, the unification of weak and elec- 
tromagnetic interactions into a renormalizable theory requires to use the principle of local 
gauge symmetry in combination with the phenomenon of spontaneous symmetry breaking 
(SSB). It is indeed the only known way to generate all the particle masses by preserving 
the underlying gauge symmetry. In the SM, one must introduce a fundamental complex 
doublet of scalar fields. However, in order to simplify the discussion, let us just consider a 
field theory with a real single scalar field (/J, as this does not alter at all the nature of the 
problem under discussion. To trigger SSB, one must introduce a potential for the field 
which in renormalizable QFT takes the form (the tree- level Higgs potential): 

y((^) = ^mV' + ^A<^^ (A>0). (2.1) 

Since we are dealing with a problem related with the CC, we must inexcusably consider 
the infiuence of gravity. To this effect, we shall conduct our investigation of the CC 
problem within the semiclassical context, i.e. from the point of view of quantum field theory 
(QFT) in curved space-time [38] . It means that we address the CC problem in a framework 
where gravity is an external gravitational field and we quantize matter fields only [14, 39]. 
The potential in equation (2.1) is given at the moment only at the classical level, but 
it will eventually acquire quantum effects generated by the matter fields themselves. In 
this context, we need to study what impact the presence of such potential may have on 
Einstein's equations both at the classical and at the quantum level. 

Einstein's field equations for the classical metric in vacuo are derived from the Einstein- 
Hilbert (EH) action with a cosmological term A^ac (hereafter the CC vacuum term). The 
EH action in vacuo reads: 

SeH = J d^X^M {R-2 Avac ) = J d^X^\ ( R - PAvac) . (2.2) 

(See the Appendix A for our notation and sign conventions). Here we have defined pAvnc, 
the energy density associated to the CC vacuum term: 

PAvac = e ''^ ■ (2.3) 
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The classical action including the scalar field ip with its potential (2.1) is 



S = Seh + / d'^x 



(2.4) 



Due to the usual interpretation of Einstein's equations as an equality between geometry 
and a matter-energy source, it is convenient to place the /OAvac term as a part of the matter 
action, S[ip]. Then the total action (2.4) can be reorganized as 



S= TTT^ / d^Xy^\R + S[^], 

lovrGiv J 



with 



S[^] 



\9\ 



1 



d^X 



(2.5) 



where C^p is the matter Lagrangian for ip. For the moment, we will treat the matter fields 
contained in as classical fields, and in particular the potential V is supposed to take the 
classical form (2.1) with no quantum corrections. If we compute the energy- momentum 
tensor of the scalar field in the presence of the vacuum term pAva.c, let us call it T^, we 
obtain 



2 5 Slip] 



where we have used d\/\g\/ dg°''^ = — (l/2)y^[^ (7afe- Here 



dab ^ip — 9ab PA vac + T"^, , 



(2.7) 



-^ab 



daP dbp - }- gab dcP d^ip 



+ gab v{p) 



(2.^ 



is the ordinary energy-momentum tensor of the scalar field ip. 

In the vacuum (i.e. in the ground state of p) there is no kinetic energy, so that the first 
term on the r.h.s of (2.8) does not contribute in that state. Only the potential may take 
a non-vanishing vacuum expectation value, which we may call {V{ip)). Thus, the ground 
state value of (2.7) is 



abi 



5abyOAvac + {T^^) = 9ab (PAvac + {V{p))) = p 



.cl 

vac yab ■ 



(2.9) 



where p'^^^^ is the classical vacuum energy in the presence of the field 93. 



If > in equation (2.1), then (93) = =^ 
energy is just the original pAvac term, 

C^ab) = 9ab PAv 



{V{p)) = and the classical vacuum 



(2.10) 



This result also applies in the free field theory case. However, if the phenomenon of SSB is 
active, which precisely occurs when m? < 0, we have a non-trivial ground-state value for 
(p: 



-6 
A 



(2.11) 
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In this case, there is an induced part of the vacuum energy at the classical level owing to 
the electroweak phase transition generated by the Higgs potential. This transition induces 
a non- vanishing contribution to the cosmological term which is usually called the "induced 
CC". At the classical level, it is given by 

PAind = {V{^)) = = {^f = -A_MlMl, (2.12) 

In the last two equalities, we have used the physical Higgs mass squared = —2m? > 0. 
Indeed, if we redefine the Higgs field as % = (p — {f), then its value at the minimum will 
obviously be zero. This is the standard position for the ground state of the field before 
doing perturbation theory. The physical mass is just determined by the oscillations of 71 
around this minimum, i.e. it follows from the second derivative of ^ at = ((/?). We have 

— 1/2 

also introduced the so-called Fermi's scale Mp = Gp ~ 290 GeV, which is defined from 
Fermi's constant obtained from muon decay, Gp ~ 1.166 x 10^^ GeV^^. 

In view of the previous SSB contribution, it is clear that we must replace — ?■ in 
the expression of Einstein's equations in vacuo. Furthermore, in the presence of incoherent 
matter contributions (e.g. from dust and radiation) described by a perfect fluid we have 
the additional term Tab = {p + p)uaUb —pQab- Therefore, Einstein's equations in terms 
of coherent and incoherent contributions of matter, plus the vacuum energy of the fields, 
finally read 

Rab - ^gabR = -8vr Gn [{Tab) + ^a^) = "StT Gn [Oab (/^Avac + PAind) + Tab] ■ (2.13) 

We conclude that the "physical value" of the CC, at this stage, is not just the original term 
PA vac, but 

PAph = PAvac + PAind , (2-14) 

where the induced part is given by (2.12). However, it is pretty obvious that a severe fine 
tuning problem is conjured in equation (2.14) when we compare theory and experiment. 
Indeed, the lowest order contribution from the Higgs potential, as given by equation (2.12), 
is already much larger than the observational value of the CC. Using the LEP lower bound 
on the Higgs mass (M^ > 114 GeV), equation (2.12) yields pAind — — lO^GeV^. Roughly 
speaking, the VEV of the Higgs potential is (naturally) in the ballpark of the fourth power of 
the electroweak VEV, i.e. pAmd ~ v'^, where v = 0(100) GeV. Thus, being the CC observed 
value of order p\ ~ 10~^^ GeV^, the electroweak vacuum energy density is predicted to be 
55 orders of magnitude larger than p^! 

Suppose that the induced result would exactly be pAmd = —10^ GeV^ and that the 
vacuum density would exactly be p\ = +10"'^''' GeV^. In such case one would have to 
choose the vacuum term pAvac in equation (2.14) with a precision of 55 decimal places in 
order to fulfill the equation 

10-^^ GeV^ = PAvac + PAind = PAvac " 10^ GcV^ . (2.15) 

This is of course the famous fine-tuning problem. This problem is in no way privative of the 
cosmological constant approach to the DE, but it is virtually present in any known model of 
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the DE, in particular also in the quintessence approach [13]. Indeed, the quintessence scalar 
field potential V{ip) is supposed to precisely match the value of the measured DE density 
at present starting from a high energy scale, usually some GUT scale (p = Mx ~ 10^^ 
GeV. In order to achieve this, an ugly fine-tuning of its initial value {V) ~ 10^^ GeV^ is 
unavoidable. Therefore, the quintessence approaches, apart from introducing extremely 
unnatural small mass scales of the order of the Hubble function at present (hence masses 
as small as m^, ~ 10"'^^ eV), are plagued with fine-tuning problems in no lesser degree 
than the original CC problem itself. 

However, this is not quite the end of the story yet. In QFT the induced value of 
the vacuum energy is much more complicated than just the simple result (2.12), and the 
fine-tuning problem is much more cumbersome than the one expressed in equation (2.15), 
see the Appendix B for a more detailed exposition. At the end of the day, we really 
need some mechanism that is able to concoct the tuning "dynamically", i.e. automatically, 
and without requiring the intervention of some carefully "designed" counterterm encoding 
the aforementioned fabulous numerical precision. In the rest of the paper, we will try to 
convince the reader that such mechanism to avoid fine-tuning does exist. 

3. J-'{R, ^)-cosmology 

As already mentioned, in Ref. [28] we investigated a powerful mechanism for relaxing 
dynamically the vacuum energy or effective cosmological constant. For all of its virtues, it is 
nevertheless based on introducing a direct modification of the gravitational part at the level 
of the field equations, and hence without any obvious connection with an action principle. 
A deeper step in the theoretical construction process would be to implement it from an 
action functional. For this reason, we investigate here the CC relaxation mechanism within 
the context of the J-{R, Q) modified gravity setup and in the metric formalism, where J- 
defines a functional of the Ricci scalar R and the GauB-Bonnet invariant 



(For more details on notation, see once more the Appendix A). The reason why the higher 
curvature term is Q rather than the individual higher derivative components that define it 
will become clear later. 

3.1 The generic class of the AXCDM models 

Our theoretical J-{R,Q) construct is placed within the general class of the so-called AXCDM 
models of the cosmological evolution [24, 28], in which the cosmological term is supple- 
mented [not replaced!) with an effective entity X at the level of the field equations. The 
resulting cosmological system is thus characterized by a compound dark energy made out 
of a (constant or variable) cosmological term p\ and the contribution px corresponding to 
a new entity X (called the "cosmon"). The total DE density (or "effective CC density") 
reads 



G = R^- mabR"^ + RahcdR 



>abcd 



(3.1) 



PAcff = PA + PX ■ 



(3.2) 
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In the simplest formulation of the AXCDM model, in which Newton's coupling Gn is 
constant, this overall DE density is locally and covariantly conserved [24]: 



PAcff + 3i7 (pAcff +PAcff) = 0, (3.3) 

where H{t) = d/a is the expansion rate (overdots represent time-derivatives with respect 
to the cosmic time). We can reexpress this equation in terms of cjefr, the effective EOS of 
the compound DE system: 

PAcS + 3 H [1 + LOcs) PAcS = , Wcff = = , , (3.4) 

PAcS PA + PX 

where = —pA is the CC pressure and px is the pressure of the cosmon component. It 
follows that matter is also covariantly conserved. 

As emphasized in [24], the cosmon is not to be viewed in general as a field, but as 
a complex object emerging from the full structure of the effective action. Its dynamics 
is determined either by suggesting an explicit modification of the effective action, or by 
providing a particular evolution law for the CC term, pA = PA{i)- In either case px = Px{t) 
is then completely determined by the local conservation law (3.3). This can be further 
appreciated by rewriting that law as follows: 

PA + PX+ 3//(l + a;x)px = 0, (3.5) 

where we have defined the effective EOS of the cosmon: ux = Px/ Px- The effective EOS 
of the compound DE system can now be written as 

-PA + ^x px 1 , /I , \ a a\ 

Wcff = , = -l + (l+a;x) . (3.6) 

PA + PX PAcff 

In general, both ujx and Wcff will be non-trivial functions of the cosmic time or the cosmo- 
logical scale factor, a = a{t), or the redshift: z = (1 — a) /a. In the original AXCDM model 
of Ref. [24], one starts from a given evolution law p\ = pa{H) for the CC term, specifically 
one which is motivated by the presence of quantum corrections, and the EOS parame- 
ter LVx for the cosmon is taken constant. This setup enabled a full analytical treatment 
and an explicit determination of the cosmon energy density as a function of the redshift, 
Px = Px{z)- In the present case, however, pA is just a constant (given by the initial p^), 
and the dynamics of the cosmon - represented by the conservation law (3.5) - will be com- 
pletely controlled by the influence of a new gravitational term that modifies the EH action. 
In fact, the complete action functional is composed, apart from the conventional matter 
part and the EH term, also of the aforementioned function of the invariants R and Q - 
whose contribution to the total action will, for brevity sake, be referred to as the ^^J-{R, Q)- 
functional". This functional determines the dynamics of X. Notice that the apparent 
simplification produced by the fact that pA is now just a constant is compensated by the 
fact that the cosmon EOS will be a complicated function of time or redshift: ojx = wx(-z)- 
As a result, a complete analytical treatment will not be possible. The reward, however, 
will be a powerful formulation of the AXCDM model in which the cosmon will be able to 
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efficiently deal (as its name intends to suggest) with the old CC problem, and specifically 
with the toughest aspect of it: the fine-tuning problem. 

Let us note that the definition of the cosmon X in the wide class of AXCDM models 
is really very general because it introduces the energy density px not at the core level of 
the action, but at the level of the field equations. In the present case, our AXCDM model 
is constructed ah initio at the level of an action functional containing the piece J-{R,Q). 
This is of course a great advantage from the theoretical point of view. However, in order 
to recognize what is px in the present case, and finally unveil the corresponding /OAeg (i.e. 
the "effective CC") that results from the presence of the J-"-functional, one has to account 
for the field equations of the complete effective action. We do this in the next section. 

3.2 Action principle and field equations 

The complete effective action of our cosmological model reads 

1 



d X \/\g\ 



R- p\-F{R,g) + C, 



(3.7) 



where Cff) stands for the Lagrangian of the matter fields. Clearly, our functional constitutes 
an extension of the EH action with CC in which we have added the J-'{R,Q) part, i.e. a 
generalization of equations (2.5)-(2.6). We take within the class of functions of the form 

J^iR,g) = f3F{R,g) + A{R), (3.8) 

in which /? is a parameter, F[R,Q) a non-polynomial function (see below) and A{R) = 
02 -R^ + a-^R^... is a (low order) polynomial of R. The latter has neither linear nor indepen- 
dent term, the reason being that these terms can already be included as part of the EH 
action with CC. Notice that we do not use Q in the structure of the polynomial A because 
the first term would just be which is a topological invariant in four dimensions. As for 
F{R,Q), we take it as a function of negative mass dimension. We will usually call it the 
"i^'-term" . The canonical implementation of it is a rational function of R and Q vanishing 
as R and G are sufficiently large, say 

with B{R, G) a polynomial in R and Q. Within this canonical ansatz (which we will adopt 
for most of our considerations), the complete functional (3.8) reads 

^<^' ^' = Bile) ^ 

Since, in the FLRW metric, F{R,G) defines a functional ^{H) of the Hubble function or 
expansion rate H{t) = a/a and its time-derivatives, we can express the required condition 
as 

F{R,g) = ^H)^0 for H>H;^, (3.11) 

where H*^^ is the Hubble rate at the time when the very early radiation period sets in, hence 
after both inflation and reheating have already taken place. We shall explain the motivation 
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for this condition later on. Since F{R,Q) dies off with curvature, the parameter /3 in 
equation (3.8) must have a positive dimension of energy, namely 

P = M^ , with Af>0, (3.12) 

where = + 4 is an integer in which is equal to the mass dimension of the B{R, Q) 
polynomial, and M represents some characteristic cosmic mass scale associated to our J^- 
functional. For instance, if B{R^ Q) is a polynomial involving just quadratic terms in R 
and linear in we would have N = 8. Notice that if the polynomial A{R) reduces to the 
monomial a2R^, then 02 is a dimensionless coefficient. In this case, the only dimensionful 
scales carried by J- are f3 and those that might involve B{R, Q) in the form of monomials 
of R (respect. G) of order 3 (respect. 2) or above. The mass scale Ai should have some 
physical significance, and therefore it will be interesting to check which are the typical 
values allowed for M (depending on the choice of F) in order to implement realistically the 
relaxation mechanism in our Universe. Generalizations of = 1 /B offer more possibilities 
for the mass scale M. and will be discussed later on in Sec. 5.5. 

Another essential ingredient of the effective action (3.7) is of course the "initial" cos- 
mological constant term 

A' = 8ttGnPa, (3.13) 

in which stands for all possible vacuum energy density contributions (of arbitrary size) 
pertinent to "initial" phase transitions in the early universe, e.g. the GUT phase transitions, 
the electroweak transition, the QCD quark-gluon transition, and in general all vacuum 
energy density contributions associated to the matter fields (bosonic or fermionic) of the 
Lagrangian C^. For instance, p\ embodies the important electroweak contribution from 
the SM Higgs potential, i.e. |pew| ~ Mp My^ ~ lO^GeV^. In addition, p\ integrates the 
QCD contribution, which is of order /Cqcd ~ ^qcd^ with Aqcd ~ 0.1 GeV. Taken alone, 
any of these SM contributions from particle physics is enormous as compared to the current 
value of the CC density, ~ 10"^'' GeV^. On the other hand, p\ may contain much larger 
contributions; in general it will be dominated by the maximum effect prevailing in the early 
universe, which should be determined by the strongest GUT phase transition, say p\ ~ M^, 
with Mx ~ 10^^ GeV. At the same time, p\ includes the vacuum term pAvaci too, if only 
for renormalizability reasons. This term is completely free and, in the traditional approach 
explained in section 2 and Appendix B, it can be used (after renormalization) to fine-tune 
all the other contributions. 

For all of its non-trivial composition, is not yet the physical CC of the J-{R, G)- 
cosmology, even though it contains all the ingredients of the traditional approach, i.e. all 
the terms on the r.h.s. of equation (2.14). If it were, the tuning that we ought to apply to 
PAvac to compensate for the energy released during the GUT phase transition would make 
the 55-digits-electroweak-tuning described in section 2 pale in comparison, and hence would 
further increase the severity of the CC problem to the utmost level! But this is not what we 
shall assume here, so pAvac will be taken as any other contribution, it is not important which 
one. What is important is that we will not take a special (fine-tuned) value for it (in contrast 
to the traditional approach, see Appendix B, and we will need not do it at any stage. 
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How to cut off from tlie root the unending escalade of fine-tunings plaguing the 
traditional approach and still render a sound value for the physical CC? A first key ap- 
pears when we compute the "effective Einstein's equations" of the J-{R, ^)-cosmology. 
They emerge from functionally differentiating (3.7) with respect to the metric, with the 
result (cf. Appendix A for details): 







5S 



6g- 



ab 



Gab + 9ab Pa + 2-E'a6 + ^aft 



(3.14) 



where the surface terms have been omitted. Therefore, instead of (2.13) we now have 



(3.15) 



Apart from the Einstein tensor G"^ = R""^ — {1/2) 6\R, we have the constant vacuum 
energy and the energy-momentum tensor Tab of matter. Notice that we have already 
absorbed in the vacuum effects (e.g. phase transitions and quantum effects) associated 
to the coherent matter contributions from the fields (j) in the Lagrangian. Therefore, Tab in 
equation (3.14) involves only the incoherent matter contributions. Finally, there is a new 
("extra") gravitational tensor E""^ coming solely from the J-"(-R, ^)-term in the action (3.7), 
which we have also placed on the r.h.s. of the field equations in order not to distort the 
standard Einstein part, usually positioned on the l.h.s. 

On a spatially flat FLRW background with line element ds^ = dt^ — a?'{t)dx^ , scale fac- 
tor a{t) and expansion rate H{t) = d/a, the tensor components of the various gravitational 
parts in (3.15) are given by 



(3.16) 



and 



E' 



E\ 



- \f{R, g) - Q{H + H'^)J='^ + GHT^ - 24H^{H + H^)J='^ + 24H^f^] (3.17) 



T{R, G) - 2{H + 3H'^)T^ + 4HT^ + 2T^ 



(3.18) 



where J-"^ = dF/dY are the partial derivatives of J- with respect to y = R,Q. Let 
us remark that in all these expressions the time derivative of the expansion rate can be 
reexpressed as H = —H^(q + 1), where 



aa 

'72" 



(3.19) 



is the deceleration parameter. This quantity will be very important in our discussions, as 
we shall see immediately. For one thing the two fundamental curvature invariants on which 
our action functional depends can just be expressed in the FLRW metric in terms of H 
and q as follows: 



R 



-2AH^q . 



(3.20) 
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Notice from equations (3.15) and (3.17) that the effective CC density in our J^{R,0)- 
cosmology, and therefore the quantity playing the role of DE in our framework, is not just 
the parameter p^, but the full expression 

PA,s{H)= Pi + PindiH), (3.21) 

in which 

And = 2£;°o (3-22) 

constitutes that part of the effective CC which is genuinely induced by the J^-functional. 
The expression (3.22) can be called the (gravitationally) "induced dark energy" ^. It adds 
up to the original cosmological constant p\ to produce the quantity PAcff or "effective 
vacuum energy density". The induced DE is obviously dynamical, and with it the total 
effective DE density too. Therefore, p\cS defined by (3.21) runs with the expansion of the 
universe. Because of (3.20), PAeff is a function of the expansion rate H, the deceleration 
parameter q and its first time derivative: pAeS = PAcs{H, q, q), but for simplicity we shall 
indicate it sometimes simply as pAcfr = PAes{H) - as we shall do also with other cosmolog- 
ical quantities. Most important, in the context of the J-^{R, ^)-cosmology, the sum (3.21) is 
the very observable quantity that should be accessible to observation, as it is this quantity 
that clearly takes the role of the CC in the effective Einstein's equations emerging from the 
action functional (3.7). Put another way, pAeff is the truly "observable DE density" of the 
J-"(i?, ^)-cosmology. Notice that we cannot disentangle observationally the two terms in the 
sum (3.21), and therefore it does not matter if is very large provided the induced term 
Pind is also large, but with opposite sign, such that the sum (3.21) leaves a small remainder. 
Obviously, for this cosmology to be realistic, we expect that this remainder is small enough 
and moreover runs only mildly with H, such that it can mimic approximately the ACDM 
concordance model. But, at the same time, and in order to avoid the fine-tuning problem, 
there must be a non-trivial dynamical interplay between the two terms in (3.21), leav- 
ing just a mild running residue at all times of the cosmological evolution after inflation, 
not just now. Let us remark that models with mildly running cosmological term provide 
a global fit to LSS and CMB data perfectly comparable to the ACDM model [22]. 

3.3 Searching for the class of T functionals 

We will argue that for a realistic approach of the new J-{R, ^)-cosmology, we expect the 
following two conditions to occur: 

i) The effective cosmological term pAcff = PAcs{H) in equation (3.21) must essentially 
coincide with the enormous value of the vacuum energy density in the vicinity of 
the post-inflationary time, when the seeds of a huge permanent vacuum energy are 

^The name seems appropriate as long as the gravitational functional induced by J-{R, Q) is treated as a 
part of the full energy momentum tensor. It would be misleading to call pind just "induced CC", because as 
we have seen in section 2 (cf. also Appendix B) this name is usually reserved for the classical and quantum 
contributions to the CC emerging from the matter part (e.g. the quantum corrected VEV of the Higgs 
potential), which we have already absorbed in p\ from the very beginning. 



-14- 



first sowed - and need to be removed. This value is wliat defines our big "initial" 
cosmological constant and corresponds to an epoch characterized by an expansion 
rate H < Mx- The function (3.8) must, therefore, allow for the following behavior 
of the quantity (3.21): 

PAcff(^) ^P\< M% {H;,^ <H<Mx). (3.23) 

That this relation can be amply fulfilled in our framework, can be argued as follows. 
To start with, remember that we have imposed the condition (3.11) on the F-term. 
For H ~ Mx such condition should presumably be satisfied (because H*^^ <^ Mx), 
and hence in this range we have 

J^{R, g) A{R) = a2R^ + ... (3.24) 

where we have taken the simplest non-trivial possibility for the polynomial A{R) in 
(3.8). Despite that this term grows with R, its value near the startup of the radiation 
epoch is not sufficiently large yet as to distort the goodness of the condition (3.23). 
For instance, for a typical GUT phase transition with A'lx ~ 10^^ GeV we have 
R ^ ~ 87r px/Mp ~ Svr M^/Mj,, and thus 



-64^' ^-10-', (3.25) 



Pa M'p 

where Mp = G~^''^ ~ lO^^ GeV is the Planck mass. Clearly, the condition (3.25) 
would not be so easily satisfied if the vacuum energy left after primordial inflation 
would be very close to Mp. However, it does not seem realistic (not even necessary) to 
try to extrapolate cosmology up to this point, so we avoid this speculative situation 
which would probably require a deeper knowledge of the space-time structure at 
the level of Quantum Gravity rather than QFT in curved space-time (as we are 
dealing with in our approach). In short, the relaxation mechanism should efficiently 
wash out the large contributions to the vacuum energy only after inflation has fully 
accomplished its role, and more specifically after the reheating processes have been 
able to "restore" the "initial" relativistic matter content of the universe in the form 
of the so-called "radiation epoch" ; but of course not before, since otherwise inflation 
itself could not have occurred. It follows that for the practical study of the relaxation 
mechanism we can simplify the induced DE (3.22) to the reduced form 



Pind PF = Pind 



(for H <Mx), (3.26) 

A{R)=0 



where the notation pp reminds us that this part is totally attributed to the F-term 
with no contributions from the polynomial A(R) in equation (3.8). 

ii) At the same time, the effective quantity p\ef[{H) must not disturb the standard 
thermal history of the universe, and should reach the present epoch with a value 
PAesiH = Hq) very close to p^ ~ 10~^'''GeV^. In view of the fact that p^ <^ 
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it means that we need a huge dynamical cancelation between the two terms on the 
r.h.s. of equation (3.21) during both the radiation and matter epochs: pind — ~Pa- 
On the face of (3.26), in practice this means that we must have 

PF ~ -p\ (y H < H*^^ < Mx until the present) . (3.27) 

This is of course the most delicate point of our construction and relies significantly 
on a suitable choice of the i^-term. For example, the choice in equation (3.9) is 
convenient because it makes allowance for the requirement 

B{R, g)^0 H < H*^^ < Mx until the present) , (3.28) 

as a starting point to fulfill the relation (3.27). Indeed, let us note that the condition 
(3.28) insures that F{R,Q) and its derivatives F^^ (Y = R,g) become arbitrarily 
large; in fact as large as p\ might be, but not infinite because p\ is anyway finite and 
hence the point i? = is actually never reached. How to make the requirement (3.28) 
natural (without fine-tuning) such that the relation (3.27) is fulfilled for arbitrarily 
large p\, is something that we will discuss in much of the remainder of this paper. 

Some further comments are now in order. The aforementioned conditions are interesting 
in that they not only define the range of the history of the early universe in which the 
relaxation mechanism of the vacuum energy must operate, they also make clear that the 
inflationary scenario can be preserved. This setup might actually prepare the ground for 
triggering primordial inflation itself through an i?^-term (or a higher order polynomial) in 
the functional A(R) in equation (3.10). Indeed, equation (3.24) is fulfllled for H > Mx, so 
that R^ becomes dominant in the far UV regime. This is consistent with the fact that the 
renormalizable quantum theory of matter flelds on a curved background must necessarily 
include the action of vacuum, which contains the higher order i?^-curvature terms [38]. 
While these terms are irrelevant for scales of order Mx or below - cf. equation (3.25) - 
they nevertheless become dominant near the Planck scale, where R ~ Mp, and in fact 
they then furnish the driving force for i?^-inflation. In other words, the J^(-R, ^)-cosmology 
provides also a possible natural connection with Starobinsky's inflation [40] and the more 
recent developments on anomaly- induced inflation [41, 42, 43]. 

Finally, a possible additional bonus of the polynomial term in (3.8) is that it may 
provide an escape to some instability issues discovered in extended gravity theories in 
the metric formalism [36] , as it is known e.g. for T{R) and J^{R, Rab^ -^abcd) theories [44] . 
There are also potential difliculties as far as concerns the astrophysical implications on 
the solar system measurements [36] . In our case, also the functional defined by J-{R,Q) 
may not be completely free from typical modified gravity problems [45]. Nevertheless, 
it behaves better than the general J^(i?, iZ^^, iZ^^^^)-functionals and this should suffice for 
the illustration purposes of our paper. For example, a remarkable feature of the J-{R,Q) 
functionals is that they avoid the Ostrogradski instability, i.e. they do not involve vacuum 
states of negative energy, and are thought to have a reasonable behavior in the solar system 
limit [46, 47]. At the same time, Gaui3-Bonnet models are automatically free from graviton 
ghosts and other singularities [47, 48] . 
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From the foregoing considerations, it became clear that the polynomial contribution 
A{R) is relevant for restoring renormalizability at high energies, for bridging the relaxation 
framework with inflation and maybe also to help curing stability issues, but the polynomial 
term is not indispensable for the relaxation mechanism itself. Therefore, we shall hereafter 
substitute the full J-{R, Q) functional (3.8) for its reduced part (3 F{R, Q), where the F-term 
is given by (3.9). 

3.4 Effective equation of state 

The matter sector can be adequately described by a perfect fluid stress tensor Tab = 
(p + p)uaUb — PQab with four- Velocity Ua-, energy density p and (isotropic) pressure p, 
respectively. For a fluid at rest [u'^ = 5q) this means 

rO=p, T) = -p5\. (3.29) 

Using a fluid description also for the vacuum energy, we define the (effective) pressures 
corresponding to the induced part (3.22) through pp = —{Tf)1/3, where from (3.15) we 
have {TfY- = 2E'j. Thus 

2 

PF = -^E\, PAcS=Pf-Pa^ (3-30) 

where the last expression is the corresponding pressure for the full effective vacuum energy 
(3.21). Since the Einstein tensor Gab and the extra gravitational tensor Eab are both co- 
variantly conserved, it follows that the J-{R, ^)-cosmology automatically conserves matter, 
too, i.e. V^Tab = 0. Therefore, the Bianchi identity on the FLRW background leads to the 
local covariant conservation laws 

Prr + 3H{pn+Pn)=0, (3.31) 

which are valid for all the individual components (n = radiation, matter and DE) with en- 
ergy density pn and pressure pn- In particular, the conserved matter {pm = 0) and radiation 
(Pr = /5r/3) components integrate immediately and yield the usual expressions pm = P%,o.~^ 
and pr = PrO-'"^, whereas the effective CC density is related to the corresponding pressure 
through equation (3.3), which in the present case is equivalent to 

PF + 3H{pF + Pf) = PF + 3H {1 + ojFit)) PF = 0, (3.32) 

because the constant part of the effective CC satisfies p\ = —p\ and therefore cancels out 
in (3.32). In the previous equation, we have defined 

^F{t) = ^ , (3.33) 

PF[t) 

which plays the role of the EOS for the F-term. As it is obvious, the conservation law 
(3.32) is a particular case of the general AXCDM one (3.5) for PA = p\ = const. The role 
of the cosmon density px is thus played by pF- The local conservation law (3.32) cannot 
be directly integrated because the EOS parameter (3.33) is actually a non-trivial function 
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of the cosmological evolution. The parallehsm with the AXCDM model can be made even 
more manifest if we define the overall EOS of the compound DE system formed by the 
constant p\ and the induced pp component as PAeff = '^eflf PAeff- It is easy to see that ujes 
is related with (3.33) as follows (for convenience we trade time for redshift): 

u,s{z) = -1 + {1 + uf{z)) //^^^ (3.34) 

Pa + PF[,z) 

which is again a particular case of the general AXCDM effective EOS defined in equation 
(3.6), with the correspondence px Pf-, — ^ and p\ + prit) = PAcS- 

The expressions for pp and pp computed in section 3.2 are not just opposite in sign. 
Thus, we do not expect LOp = —1, and therefore Wcfr 7^ ~1 either. Both are non-trivial 
functions of time or redshift: up = iOp{z) and Weff = u}es{z). Models with variable cosmo- 
logical parameters indeed usually exhibit this feature [49]. In our model universe, we see 
from (3.22) and (3.30) that the departure from a strict cosmological constant behavior is 
caused by the fact that E^q ^ {1/3)E^^ in equations (3.17) and (3.18). This relation would 
hold only if all derivatives = dF/dY = {Y = R,Q), but this is impossible for the 
typically needed structure for F, see equation (3.9). The uJcs — 1 feature, therefore, 
will actually persist for the entire cosmological history, but we expect that the departure 
from —1 will not be very important near our time because we are currently observing a 
predominance of the "DE epoch" (p^, pr <^ PAcs)i in which the DE behaves essentially as 
a CC term. So we should ensure that a;efr(-z) — —1 for redshift z ~ 0. 

Well within the original spirit of the class of AXCDM models, the local covariant con- 
servation law satisfied by the cosmon - equation (3.32) - plays a fundamental role to 
elucidate the dynamics of the model. Thanks to this conservation law (which acts as a 
first integral of our dynamical system), there is no need to use the complicated expression 
(3.18) - which leads to a differential equation of one order higher than (3.17). Therefore, 
it is not necessary to use (3.18) to find pAcff through (3.30). In practice, the effective EOS 
can be determined with the help of the local covariant conservation law (3.3) or just (3.32). 
For this, PAeff is to be determined first. 

But how to find explicitly the effective vacuum energy pAcff? To this end let us consider 
the "Q-component of the Einstein equation (3.15): 

SH"^ = SirGNiPm + Pr + PAcff) • (3-35) 

This is the generalized Friedmann's equation for the -F(i?, ^)-cosmology, in which the role 
of the CC is played by pAcff- Such equation provides the clue for integrating the field 
equations and determine all relevant energy densities. Substituting equation (3.22) in 
the formula (3.21) for the effective vacuum energy pAeff, one obtains an expression that 
depends on H, q and q. In practice, for any function /, we can trade the time evolution 
of the resulting expression for the scale factor dependence through / = aH[a) f'{a). In 
particular, H{a) = aH{a)H' (a) and moreover from q{a) = —1 — aH'{a)/H[a) we have 
q{a) = —aH'{a) — a?[H"{a) — H''^{a)/ H{a)]. Proceeding in this way with all the terms, we 
arrive at the form pAcff = PAes{cL, H{a), H' (a), H" {a)). Therefore, Friedmann's equation 
(3.35) can be finally cast as a second order differential equation for H{a) that can be 
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solved numerically. Once H{a) is known, the effective vacuum energy p\eS = PAesia) is also 
known, and can be plugged in the local conservation law (3.3) to determine pAcS = PAefr(a), 
and from here the desired EOS Weff(a) = PAcfr('3^)//'Aeff(a) ensues. We shall follow this 
procedure in practice. 

Furthermore, the '^-component contains information on the acceleration, and can be 
expressed in terms of the deceleration parameter (3.19) as follows: 

SH"^ q = AttGn [2 pr + Pm + (1 + 3a;efr) pacAf] • (3.36) 

Using (3.35), we can recast this expression such that its r.h.s. contains just the sum of 
pressure components: 

H^{q - 1/2) = 47rGN{Pr + PAcff) = 47rGAr [pr/3 + a;cfr pAcff] • (3.37) 

Similarly, we obtain 

3 H'^{1 - q) = AttGn [pm + PAcff - 3pAcff] = AttGn {pm + (1 - 3u;cff) PAeff) ■ (3.38) 

The last two equations are convenient forms of the dynamical equation for the acceleration, 
they are written in terms of the deceleration parameter (3.19) and will be used later on. 
By ignoring pAcff for the moment it is easy to see from (3.38) that the radiation epoch (in 
which Pm can be neglected) is characterized by g = 1, whereas from (3.37) it is transparent 
that the matter epoch (in which p,. can be neglected) is characterized hy q = 1/2. In 
Sees. 5.1 and 5.2 we will show that these observations still hold when taking into account 

PAeff- 

Let us note that the above formulae just follow the normal pattern of equations char- 
acterizing a cosmological medium which is composed of several fluids. If these fluids have 
EOS parameters = Pn/Pn and density parameters ^n{o) = Pn{o)/Pc{a) normalized 
with respect to the critical density Pc(o) = 3H'^{a)/{8'irGN), one can easily show that the 
deceleration parameter (3.19) can be expressed as 

g = (1 + 3uJn) ^ ; Y.nn = l. (3.39) 

n n 

Clearly q = (1, 1/2, —1) for radiation {ljr = 1/3), matter (cjm = 0) and standard vacuum 
energy (cja = ~1) dominated epochs respectively, which correspond to having the dominant 
density parameter in each epoch 0^ = 1 and all the others zero. In the present epoch, we 
have a mixture of matter and DE in which the latter behaves very approximately as vacuum 
energy, therefore the current value of the deceleration parameter in the ACDM model is 
generally expressed as qo = r2^/2 — O^. The only note of caution is that, within the 
framework under consideration, the effective CC term pAcff does not behave as standard 
vacuum energy because it has a non-trivial EOS lOcS, which is not equal to —1, and in 
general is a complicated function of time cjgff = ^cs{t) or of the redshift ~ see (3.34). 
With this only proviso, equation (3.36) is easily seen to follow from the general one (3.39) 
accounting for a mixture of fluids. In this way, defining Q^csiO') = PAcs{cl)/ Pc{ci) also for 
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the effective vacuum fluid of our model, the value of the deceleration parameter at present 
is to be written as 

qo = + (1 + 3a;°ff)OLff/2 = ^ i'^ + ^c^efF ^^Aeff) , (3-40) 

where the second equality is valid only if the universe is spatially flat. Here and 
cUgjj are the current values of these quantities. Equation (3.40) is obviously consistent with 
(3.36) when the radiation contribution is neglected. In sections 5.1 and 5.2 we will discuss 
the relation of the non-trivial EOS lOcs of the effective vacuum energy pAcfr with the EOS 
of matter and radiation in the various epochs. 

3.5 Evading a "no-go theorem" 

Before further exploring our J-{R, ^)-cosmology, let us briefly comment why it has a chance 
to evade Weinberg's "no-go theorem" for dynamical adjustment mechanisms of the cosmo- 
logical term [3]. The theorem is formulated for a system of scalar fields <pj non- minimally 
coupled to gravity, and the basic claim is that it is impossible to find a stable vacuum state 
for this system that coincides with flat space-time in the gravity sector, unless fine-tuning 
is used. In a very simplified form, the proof is based on studying the consistency of the 
combined set of equations defining the existence of the necessary extremum for constant 
scalar fields and metric, namely the set of derivatives of the matter Lagrangian with respect 
to all the fields equated to zero: 

-^ = 0, T^=0. 3.41 

O^j Ogab 

The solutions of this system must be compatible with the solution of the trace T = T"^ of 
the energy-momentum tensor being zero at the same point in (constant) field space, i.e. 

T = 2gabj^-^C^ = 0. (3.42) 

What are the chances for this possibility? A good start would be that the two expressions 
on the respective l.h.s. of the system (3.41) would be proportional, or related by a linear 
transformation. Then, a ground state solution in yjj-space would automatically be compat- 
ible with a solution of constant metric gab, which we may suggestively call rjab (Minkowski). 
Besides, this would also imply that the first term on the r.h.s. of (3.42) would vanish. Un- 
fortunately, this does not guarantee yet that the energy-momentum tensor vanishes unless 
the second term, viz. — 4^^ = +4Vcf[{^j), also vanishes. But this could not happen unless 
we would fine-tune to zero the value of the ground state of the effective potential of the 
scalar fields, quite in the same contrived way as discussed in section 2 - except that here 
we would have exactly zero on the l.h.s. of equation (2.15). Therefore, in general there is 
a "no-go" conclusion about the possibility of having r]ab as the metric solution just at the 
point of field space where it is localized the ground state of the tpj fields. In the old days, 
it was expected that if this ground state value is zero, then there would be some hope that 
some symmetry or dynamical mechanism would help reaching it without fine-tuning. But 
it does not seem to be the case, as Weinberg's no-go theorem claims [3]. 
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The root of the problem hes on the fact that the above system of equations is over- 
constrained. If we would, instead, not require to have a constant solution rjab in the current 
vacuum state, the problem should not arise, because then the metric at any time - and in 
particular in the present universe - could be a dynamical one, typically the FLRW metric, 
with gab sl function of the scale factor a = a{t). In this case, the second equation in (3.41) 
would not hold. If, in addition, we do not require that the expression for the trace of the 
energy-momentum tensor vanishes at the ground state for matter fields - and in particular 
neither at a point where gab is constant - the system becomes less and less constrained and 
we should not expect any impediment to demand that T carries, at the present time, some 
non-vanishing, even if small, energy density compatible with the curved space-time metric 
of our current epoch. In short, by loosing the constraints (by allowing dynamical metric - 
hence space-time curvature - and non-zero vacuum energy at any time) the no-go conclusion 
disappears. This is exactly our situation. The "only" final difficulty lies in achieving the 
right non-vanishing value of the current T, namely one which is small enough for particle 
physics standards. Here is precisely where the full power of the relaxation mechanism 
enters. The rest of the paper is devoted to explain why and how this is possible. 

4. Dynamical relaxation of the vacuum energy 

In Einstein's General Relativity, the theoretically expected large vacuum energy density p\ 
which was released at the early stages of the cosmic evolution would drastically change the 
essential features of the standard cosmological paradigm, in particular it would prevent the 
well-established thermal history and all the astounding successes of the Big Bang universe. 
This problem can be solved either by extreme fine-tuning or by a dynamical CC relaxation 
mechanism, which is the subject of this work. 

The big value p\ should prevail at times prior to the radiation epoch, in particular 
during the fast de Sitter expansion that characterizes the primordial inflationary phase. 
However, a "residual" vacuum energy of respectable (even of comparable) size is expected 
to remain in the universe in the vicinity of the incipient radiation epoch, i.e. the epoch 
that ensues after the universe loiters for a while in the reheating state, namely that state 
which is responsible for "re-creating" all the (relativistic) matter out of the decay of the 
inflaton or any other inflationary driving force. In fact, there is no reason to expect that 
after inflation the universe will roll down into a vacuum state of very small energy. The 
"residual" energy left in the reheating vacuum can be called again p\, because it could 
perfectly be of the same order of magnitude. Let us take into account that nothing is 
accessible to us before this time, and much less if we move deep into the inflationary era. 
Therefore, the relaxation mechanism must be operative only after inflation has ceased 
and the turbulent state of the universe, caused by the reheating mechanism, has flnally 
homogenized the fluid and triggered the primeval radiation epoch within the FLRW metric. 

Of course we cannot easily describe the interpolation processes that made possible 
the transition from the de Sitter inflationary phase into the FLRW phase, and much less 
without a fundamental microscopic understanding of the very early universe (string theory, 
brane- world, M-theory?). This goes beyond our main purpose in this paper, which is only 
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to demonstrate that a dynamical mechanism to relax the CC can be explicitly constructed. 
Thus, we shall just assume that the transition took place and that, after reheating, the 
universe was left in principle with a significant vacuum energy of the order of the initial 
de Sitter one. There is no reason whatsoever (apart from an unacceptable fine-tuning 
of the initial conditions) to expect that a sizeable vacuum energy is not there, so unless 
the universe unleashes automatically some countermeasures to reduce it fast at a minimum 
level, it may completely ruin the onset and full development of the standard thermal history 
of the Big Bang model, in particular the primordial and very successful nucleosynthesis of 
the light elements. For this reason, the neutralization process of p\ must be immediately 
put to work with utmost efficiency. 

The previous description tells us when our relaxation mechanism is supposed to start 
working. But the next (and highly non-trivial) question is: how does it work? To introduce 
the mechanism of relaxation in our modified gravity framework, let us make some ansatz 
within the class of the functionals J-{R,Q) defined by equation (3.10). Remember that we 
need to satisfy some properties described in section 3.3. For definiteness, let us choose for 
the polynomial in R and Q the expression 

B{R,g) = b2R^ + cg + bnR'', (4.1) 

in which n is some integer different from 2. This simple ansatz is already sufficient to 
explain the basic principle, which carries over to more complicated models like F = R? / 
etc that we will address briefly later on. It is obvious from the choice (4.1) that in general 
the corresponding F-term will trivially satisfy the condition (3.11). But a more difficult 
question is whether it can satisfy the dynamical neutralization condition (3.27) as well, 
which must hold for all H below H*^^ until the present epoch. Amazingly enough, achieving 
this feat is possible by an appropriate choice of the first two coefficients of the polynomial 
(4.1), whereas the third coefficient and the power n just control some smoothing properties 
of the thermal history, as we shall see below. Notice that if 62 and c are dimensionless, 
then bn is dimensionful, with dimension 4 — 2n of mass, i.e. M^"^". At the same time, it is 
clear from dimensional analysis that, for the present ansatz, the coefficient /3 in equation 
(3.12) has mass dimension N = 8. 

4.1 A toy model 

Let us now explain how the relaxation mechanism works, and let us do it by making use of 
a simpler (albeit non-trivial) toy-model example which contains already some of the main 
ingredients. Relaxation means that the observed energy density pAeff can be made much 
smaller in magnitude than the initial p\. This can be achieved by making the induced 
part pind (3.22) large in magnitude, and opposite in sign to p\, such that the two terms in 
equation (3.21) conspire to keep the sum \p\ + pi^^l ^ IPaI ^^^^ stretch of the post- 

inflationary cosmic evolution. With the current ansatz, the previous condition is realized 
when the denominator B is sufficiently small, S — ?• 0, but non-zero. 

As explained in section 3.3, for the study of the relaxation mechanism we can take the 
reduced form pp (3.26) of pmd- As a warm-up, let us consider the polynomial (4.1) for the 
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particular case where 62 = 
B = R"^ = mH^{l-q)^ 



1 and c = 6n = 0. Then, B takes the simplest possible structure 
see equation (3.20) - and we have 



PAefT(^^) = Pa + PF = pi + {3 



1 



36H^{l-qy 



+ SiH,q,q) 



(4.2) 



where £{H, q, q) represents the terms in (3.17) containing derivatives of J- (i.e., essentially 
of F). Notice that this toy-model example is similar to the one studied in Ref. [28], see 
equation (6) of the latter. The difference, however, is that here we use 1 /i?^ rather than 
1/R, and that we have the presence of the function £{H, q, q). The latter is a direct conse- 
quence of performing our analysis of the relaxation mechanism from an action functional 
rather than imposing the form of the new terms at the level of the field equations. However, 
none of these differences will change the qualitative behavior of the relaxation mechanism 
nor the fact that this setup, despite it contains the first clues to the dynamical relaxation, 
is still too simple for making it work realistically. 

The dynamical relaxation of pAefr originates from the (1 — (?) factor in the denominator 
of the expression pp in equation (4.2). The large p\ left over in the immediate post- 
inflationary period drives the deceleration parameter q to larger values until g — )• 1, which 
corresponds to radiation-like expansion. In other words, the very existence of the radiation 
period is triggered automatically by the presence of this term, which can be thought of 
as a countermeasure launched by the universe against the presence of the large "residual" 
vacuum energy at the pre-radiation era. In view of the form of (4.2), we expect that there 
will be a significant dynamical neutralization of during this epoch, for an appropriate 
sign of the parameter /?. Although q is driven dynamically to q' — )■ 1, it cannot cross 
q = 1 from below since, then, pp would dominate over p\ and stop the cosmic deceleration 
before q reaches 1. 

Let us also clarify that the function £ in (4.2) is not just a passive spectator, as it 
contributes alike to the neutralization process. The reason is that the terms of £ with 
derivatives and F^ furnish contributions to pp of the form 1/(1 — q)^ and 1/(1 — q)^, 
respectively. Thus, we end up with a general expression of the form 



pAcs{H) = p\ + PF = Pk + P 



M2 M-i 

(l-(7)2 + (l-g)3 + (l-g)4j 



(4.3) 



where the functions J\fi{H, q, q) (i = 2,3,4) do not contain the factor B. As we see, the 
additional dynamical terms on the r.h.s. of equation (4.3) stay on equal footing as far as 
their ability to neutralize the p\ term. As advanced in point ii) of section 3.3, one can 
show that the validity of the argument is general for any F(R, Q) of the form (3.9). At the 
same time, by dimensional reasons we have Ni{H, g, g) — t- as — ?• 00, and therefore the 
condition (3.11) is satisfied. By the same token, Mi{H, g, g) — )• 00 as — )• 0. 

Despite there is a tremendous cancelation in (4.3) between p\ and the "/3-terms", i.e. 
\p\ + Pp\ <^ IPaI' there is in fact no fine-tuning anywhere. The compensation is dynamical, 
and hence automatic, i.e. triggered by the evolution itself of the universe. The point is 
that p\ and pp want to drive the deceleration parameter q to different directions. Let us 
e.g. consider a dominant negative vacuum energy density < (as it would be e.g. the 
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case of the electroweak energy of the Higgs potential in the SM (see section 2). Then, 
PAeflf ~ /9\ < at the initial stage of the radiation epoch. This big negative vacuum 
energy would tend to produce a dramatic deceleration of the expansion, but at the same 
time q is fast driven to 1 until the terms 'va. PF that increase with inverse powers of (1 — g) 
become sufficiently big to compensate for /?^. Put another way, pp acts as a "dynamical 
counterterm" . Ultimately, the "fine-tuning" between and /9_p is indeed there, but it is 
not "man-made" , it is rather dictated dynamically by the universe itself! 

Worth emphasizing is the fact that the relaxation solution is dynamically stable. To see 
this, take again the case of a large and negative p\. The driving of H (by p\ < 0) to small 
values becomes compensated by the large and positive contribution of pp, which, as we have 
seen, grows as H decreases. On the other hand, any attempt of H at growing inordinately 
large would be deactivated automatically by the decreasing pp, which would make the 
term < to take over again and render H stable. In other words, p\ and pp monitor 
each other, and this feedback results in the complete stabilization of the expansion rate. 
As already mentioned, this stabilization is what impedes q ever reaching the exact value 
1. It just approaches 1 the exact amount to get the p\ term sufficiently counterbalanced. 

In this dynamical relaxation process for the CC, tiny changes of the deceleration q 
near 1 are sufficient to compensate for changes in the detailed structure of F or for large 
variations in the value of p\. In the latter case, it means that the mechanism automatically 
self-adapts to any modification of the initial conditions setting the value of p\. For example, 
it works equally well if the original vacuum energy is of the order p\ ~ (10^ GeV)^ (as in 
the SM) or if it is much larger (e.g. in a typical GUT) and very much accurate, say with 
the precise value p\ = (5.648310279 x 10^^ GeV)^ etc. This self-adapting dynamics is also 
the reason for the absence of fine-tuning in our setup. 

Note that, at the equilibrium point (in this example g ~ 1), both terms in pAeS 
are almost equal to each other apart from opposite signs (for an appropriate sign choice 
of P). Therefore, each term p\ and pp could be well approximated as a cosmological 
constant. However, their sum is not constant in general, which is the result of the implicit 
time-dependence in pp{H{t)). If mild enough, the running property of p\csiH) with the 
expansion rate can remain almost undetected to us and can perfectly simulate the ACDM 
model. Overall, two large approximate CC terms conspire to give a much smaller CC-like 
term, the observed one! We shall see explicit numerical examples in section 5. Besides, 
there is a corresponding compensation of the terms pp and —p\ in the effective vacuum 
pressure PAeS in (3.30), and we have already seen in (3.34) that the EOS of the effective 
vacuum energy density is not constant in general. 

Let us point out that the cases p\ < and p\ > are qualitatively distinct. If 
p\ is large and negative (e.g. as in the electroweak vacuum), then the driving of g to 1 
is enforced automatically by the pp term in (4.2), just to avoid that < 0. In this 
sense, the p\ < vacuum "brings forth" the radiation epoch as something inevitable after 
the primordial post-inflation period. However, in the alternative situation in which p\ is 
positive and large, in principle nothing prevents the universe from still continuing in the 
de Sitter phase, unless the vacuum energy starts to decay into radiation (e.g. by virtue 
of some particle physics processes associated to the reheating mechanism). We cannot 
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describe this decaying mechanism in our framework, but we must assume it has happened, 
and hence it should trigger the formation of "bubbles" of the q = 1 state in the vacuum. 
From here onwards the relaxation mechanism takes its turn an can automatically remove 
most of the vacuum energy from this state, thereby transforming it into a a normal heat 
bath of relativistic particles. Only after most of the vacuum energy would be neutralized, 
the evolution of the relativistic particles (radiation) in the heat bath could follow the 
pattern of the standard FLRW radiation epoch. 

In summary, in this toy-model example we have all the essential ingredients for the 
relaxation mechanism to work. The latter may not only trigger the appearance of the 
radiation epoch (specially if p\ < 0) and protects it from the devastating effects of a large 
vacuum energy remnant in the post-inflationary time; quite remarkably, it also predicts 
a very small value of the effective vacuum energy at the present time, which, amazingly 
enough, is the most sought-for "miracle" needed to solve the big cosmological constant 
problem. Indeed, in the current epoch the condition g ~ 1 has long ceased to hold and the 
compensation of p\ by pp in equation (4.2) can only occur because H in the denominator 
of pp has attained a very small value. How small is this value? The presence of the 
complicated term £ may obscure an analytic estimate here, and although we shall present 
later an exact numerical solution of a more realistic model, let us now simplify things 
momentarily by considering the effect of the F-term only ~ i.e. imagine that the if^-term in 
(4.2) is absent. This is tantamount to say that we assume A/s = = in equation (4.3). 
Then, since in the current universe we must have \p\ + pf{H)\/ p\ <^ 1, with pp ~ f3/ H^, 
it follows that the value of H that solves this equation is approximately given by 

which is well defined because, in this example, < and hence we have to choose /? > 0. 
Furthermore, from this expression it is patent that the small value of at the present time 
is just caused by the large magnitude of p\ at the early times! The values of and pAcfr 
are connected by equation (3.35), i.e. approximately by 3i/^ ~ SttGn p\cS (if we neglect 
the current matter contribution, which is anyway smaller than the observed CC). Finally, 
we can attain ~ Hq by an appropriate choice of the magnitude of the parameter /3, or 
equivalently by the mass scale Ai in equation (3.12), with N = 8. For instance, taking 
p\r^ M^, with Mx ~ 10^^ GeV, and using Hq ~ 10^^^ Q^y, we easily find M ~ 10"^ 
eV, which is in the range of light neutrino masses, i.e. a reasonable mass scale for particle 
physics standards! 

4.2 More realistic cosmological models 

As we have seen in the previous section, the simple choice 62 = 1 and c = 6^ = made for 
the coefficients of the polynomial (4.1) provides a cosmology endowed of truly remarkable 
properties. Unfortunately, that choice is too simpleminded for a realistic description of our 
universe. The reason is that while the cosmos in that scenario goes through a radiation 
epoch (g = 1) there is no possibility to drive it into a subsequent matter epoch {q = 1/2). 
Obviously, this is a fundamental shortcoming. In the following we shall try to amend this 
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difficulty and we will discuss the evolution of the effective vacuum energy in a more realistic 
CC relaxation model in which all relevant epochs are finally included. This case will be 
more complicated than the toy-model from the previous section, but the working principle 
is the same and for this reason we have explained it with some detail there. Our starting 
point is a model still based on the F-term ansatz (3.9) and with a i3-polynomial of the 
generic form (4.1). Again, for the analysis of this setup it is useful to characterize the CC 
relaxation (|/9Acfr| ^ I Pa I) with the condition i? — )• 0, although B does not strictly vanish, 
as we have explained in the toy-model example. From this condition, we will derive in 
the next section approximate analytical results for the evolution of the effective CC term 
(3.21), which we shall support with numerical simulations. 

For a realistic model of this kind, the polynomial B must have appropriate coefficients 
62 and c such that the and G terms produce a neat {q — 1/2) factor, and besides we 
need a non- vanishing coefficient 6„ to insure that the i?" term will provide a (q — 1) factor 
as well. Therefore, we introduce the polynomial 

BiR,g) = + + (y^)n ^ 24i/4(^ _ _ 2) + [6yH\l -q)]\ (4.5) 

where we have used equation (3.20). Comparing with the example above, it is easy to 
see that the second term ~ H'^^{1 — g)" will relax the effective CC in the radiation era 
((; ~ 1). In order for this term to dominate at Hubble rates H characteristic of that epoch, 
we must require n > 2 as this insures that the last term of equation (4.5) increases faster 
than at high H - i.e. faster than the first term. The latter, on the other hand, will be 
responsible for the relaxation in the matter era (<? ~ 1 /2) for lower values of H, namely for 
H < Heq, where Heq is the Hubble rate just at the transition time from radiation to matter. 
Numerically, i^eq ~ Hq, corresponding to a temperature of T ~ eV. Additionally, the 
exponent n > 2 determines the smoothness of the radiation-matter transition. Finally, the 
dimensional parameter y fixes the redshift of the transition. Obviously, it will be of order 

(4.6) 

as this is the point where the two terms on the r.h.s. of (4.5) will be of the same order. 
We thus have only two free parameters in the polynomial (4.5), which if added to the 
parameter /3 (or, equivalently, the mass scale M) in (3.12), it makes a total of three free 
parameters: 

(/3,n,y). (4.7) 

With only this small number of parameters the relaxation mechanism can be made to work 
in a pretty realistic way, as we shall demonstrate explicitly in the next sections. 

Before closing this section, the following comment is in order. For H S> H*;^ (see 
section 3.2), and specially near the de Sitter phase at ~ Mx, the deceleration param- 
eter should be forced by the mechanism of primordial inflation to stay near q = —1, and 
therefore the F-term should satisfy the condition (3.11) since the polynomial (4.5) becomes 
numerically large at high H when it is away from the region where q = 1. Notwithstand- 
ing, we must admit that we do not have at present a precise control of the interpolation 
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regime between the inflationary period and the onset of the standard FLRW cosmological 
evolution. This is actually a general problem plaguing all inflationary models. Therefore, 
we are not supposed to describe at this stage the corresponding evolution of the relaxa- 
tion J^-functional from one period to the other. In particular, the functional form of the 
F-term in the general structure of J- could change during inflation; in fact, its ultimate 
origin goes beyond the scope of this investigation. But irrespective of the details of the 
underlying fundamental theory of the J-"-functional, we expect that the condition (3.11) 
should be satisfled in order to preserve the mechanism of inflation prior to the startup 
of the FLRW cosmology. At energies close to ~ Mx the behavior of the complete T- 
functional (3.8) should not interfere with this fact, and it is thus reasonable that it takes 
the polynomial form (3.24). After all, this UV form of the effective action (if no powers 
higher than i?^ are involved) is the one that is expected for the standard renormalizable 
effective action of QFT in curved space-time [38, 39]. At the same time, it would naturally 
provide Starobinsky's type mechanism of primordial inflation [40] and modified formula- 
tions thereof [41, 42, 43]. This is a most natural expectation in a framework where the 
main job of solving the cosmological constant problem is accomplished precisely by gravity 
itself rather than by introducing extraneous scalar fields. 

5. Numerical analysis of specific relaxation scenarios 

Let us now consider the detailed numerical analysis of the J-{R^ ^)-cosmology based on a 
F = l/S-term with B the polynomial given in equation (4.5). We will consider a separate 
discussion of the matter epoch, the radiation epoch and the late time epoch. It is also 
interesting to study the future behavior, in particular to analyze if it is a pure de Sitter 
phase, as in the concordance ACDM model, or there are some significant deviations from 
it. In this section, we will perform a numerical analysis of the exact equations presented 
in section 3.2 and we will compute the precise evolution of the following basic quantities: 
deceleration parameter, the EOS parameter and the density parameters for the various 
energy densities, i.e. 



where pc{z) = 3 H^{z)/8ttGn- In all cases we will present the evolution as a function of 
the cosmological redshift z or, equivalently, with the scale factor: a = 1/(1 + z). For the 
exact numerical study of the quantities (5.1), let us recall from the discussion presented 
in section 3.4 that the generalized Friedmann's equation (3.35) provides the clue for inte- 
grating the field equations, as it determines a second order differential equation for H{a) 
which can be solved numerically, and thereby all the quantities (5.1) can be accounted for 
too^. However, in order to better understand qualitatively the meaning of the numerical 
results, we will precede our numerical analysis with an approximate analytical treatment 

^The concrete examples used in Figs. 1-8 should suffice to illustrate the working ability of the relaxation 
mechanism, even though the values of p\ in realistic GUT theories are higher than those used in our 
numerical analysis. The reason for using smaller values is simply to avoid unnecessary numerical difficulties. 
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of the behavior obtained in the various epochs. As we will see, the model under consider- 
ation faithfully reproduces the standard matter and dominated epochs, in contrast to the 
traditional modified gravity models [37], and it leads to an asymptotic evolution that may 
effectively appear either in quintessence-like, de Sitter or phantom-like mask. 

5.1 The matter era and the cosmic coincidence problem 

Let us start in the matter era where 5 ~ ^ and ~ pm oc a^^ . We are assuming that 
the matter era under consideration is not too a recent one, i.e. we suppose that the matter 
density dominates over the vacuum energy {pm > PAcs)- An exception will be discussed in 
section 5.6. To compute the EOS of the DE (i.e. of the effective CC) in this epoch we can 
obtain an analytical approximation as follows. We apply directly the relaxation condition 
i3 — > in equation (4.5), which leads to H^{q — ^) oc H^"^. Furthermore, from (3.37) 

H^q - ^) = AnGNilpr + PAes) oc ^ ^-3(n-l) _ (5 2) 

Notice that this relation becomes the ACDM scaling law pr = p^ in the radiation epoch, 
only if n = |. From (5.2) we find PacS = ~^Pr + cia~^^'^~^^ with a constant ci, which 
implies 

PAcS = C2a~^ - Pr H -^■a~^^'^~-^\ C2 = const. (5.3) 

n — 2 

as a result of solving the local covariant conservation law (3.3). Therefore, the dark energy 
EOS for the effective CC in the matter epoch can be approximated by 

WcfT(a) = = , ei -3(n-l) ' ^^■^> 

which is a non-trivial one. When n = 7/3, it actually interpolates between dust matter 
(cjeff — 0) at late times and radiation (oJcS ~^ 1/3) in the early matter era. Depending 
on the integration constant C2 a pole might occur in cocs, which can be seen in some of 
our numerical examples, see e.g. Fig. 1. The term proportional to ci could lead to an 
intermediate scaling if n > 2. For larger values of n it is not important. Notice that these 
poles in the EOS have no physical significance since all physical quantities (energy density 
and pressure) are well defined at all redshifts. The pole appears only when p\eS = 0, 
but this is no real singularity as the description in terms of the EOS parameter is not 
fundamental, it is only convenient, see e.g. an analogous situation in [24]. 

It is quite interesting to remark that the EOS analysis suggests that dark energy 
behaves like dark matter in this epoch, and we will speculate on possible applications in 
Sec. 5.6. In addition, the tracking relation p\eG oc pm refiected in (5.3) - in which the 
last two terms on its r.h.s. decay faster than the first and should thus be comparatively 
negligible ~ is quite noticeable too, and can be considered as a cornerstone for solving the 
coincidence problem in the T{R, ^)-cosmology. Recall that pm = Pm ^ ^i^d also that 
C2/Pm need not be a very small number. But even if C2/Pm — 0.1 the tracking property is 
remarkable because it shows that the DE and DM densities are following parallel evolutions 
at diff'erent levels, and one is not necessarily infinitesimal as compared to the other. If so, 
this is an enlightening clue to explain the coincidence problem. 
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Figure 1: Deceleration g, dark energy EOS Woff = PAcs/pacS, and relative energy densi- 
ties fin = Pn/Pc of dark energy pAeff (orange thick curve), dark matter pm (black dashed-dotted) 
and radiation (red dashed) as functions of redshift z in the model F = 1/B with n ~ 3, 
y^O.Tx 10-3iJ-2/3^ ^ -10*50 GeV^ r2°„ = 0.27, = lO""^, qa ~ -0.6, go = -0.5 Hq. In the 
deceleration plot the thick orange curve corresponds to the modified gravity model, and the black 
dashed-dotted curve to ACDM. 

The reason is that when the universe abandons progressively the matter epoch, i.e. 
when q starts to significantly deviate from 1/2, then the condition — )• is no longer 
implemented through q — )• 1/2, but by a significant depletion in the value of the expansion 
rate itself, which tends to a very small value H ^ H^, near Hq . This breaks the approximate 
"parallel" evolutions of pm and p\eS (equivalently of Qm and i^Aes) since the universe 
becomes DE dominated. Thus the two curves must cross (the DE one emerging from below 
because it was subdominant earlier) and then they meet at some point near our present. 
This can be seen in the numerical examples presented in Fig. 1 and Fig. 5-8, which sustain 
our claim that this could be a nice possible explanation for the cosmic coincidence problem. 
But, even more remarkably, is the fact that it emerges from the same mechanism solving 
the old CC problem! 

5.2 The radiation era 

Conventionally, the radiation era is considered as the cosmological time interval which 
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began after the phase of reheating (subsequent to the period of primordial inflation) and 
ended with the beginning of the matter era, i.e. when H takes values in the range H^q S 
H < H*^. Accordingly, the deceleration is (; ~ 1 and ^ pr oc a~^. Using again the 
relaxation condition S — )■ in (4.5), which in this case amounts to [H'^{1 — q)]^ oc H'^, we 
find 

R = 6H^{l-q) (X o^a~^ . (5.5) 



Note from (3.38) that the above equation for H^{1 — q) will exactly match the scaling 
behavior ~ of the ACDM model in the matter epoch (for negligible pa) only if n = 8/3. 
Therefore, we should choose n close to this value to obtain a smooth radiation-matter 
transition close to the standard model. Equation (3.38) leads to the relation 

3H'^{1 - q) = AirGNiPm + PAcS - ^PAcs) fx a"n, {p-m = p^a"^) , (5.6) 

and thus to pAcS = gPm. + ^PAefr + 030"^ with C3 = const. The dark energy density follows 
from inserting this expression in the local covariant conservation law (3.3) and solving it, 
with the result 



PAeff = C4a - pm + 3c3 |^-4 + " J " " > (5-7) 

where C4 = const. The corresponding EOS of the effective CC in this epoch interpolates 
between radiation (cjeS ~ 1/3) at early times and dust matter (ojcfi — 0) at later times: 



Also here, a pole in ojcs might exist depending on C4, which can be seen in Fig. 1. Again if 
n >2, the term proportional to C3 could lead to an intermediate scaling, which is absent for 
larger n. Remarkably, we find the tracking property pAefi oc p,. oc also in the radiation 
era. 

An important point to be stressed in the radiation epoch is that the aforesaid tracking 
property does not endanger the primordial nucleosynthesis process. In fact, when compar- 
ing the scaling law for radiation pr = p^a~'^ with equation (5.7), nucleosynthesis requires 
that the ratio C4/p^ is sufficiently small, typically one order of magnitude at least. We 
explicitly confirm from the numerical examples in Fig. 1 and Fig. 5-8 that this is indeed 
the case in the large z region for z > 10^, which comprises in particular the nucleosynthesis 
segment around z ~ 10^. 

The tracking feature of the DE in the matter and radiation epochs, when inspected 
in Figs. 1 and 5-8, is actually better captured if we focus on the evolution of the EOS 
parameter as a function of the cosmological redshift, ojcs = ^osiz), rather than looking 
at e.g. the corresponding behavior of the density parameters 0^. The point is that, from 
the figures one might get the wrong perception that il-AcS traces fi^ rather than ilrn in the 
matter epoch, and vice versa in the radiation epoch. This is only due to the numerical 
smallness of ^AeS in both the matter and radiation epochs. Instead, when we look at the 
EOS plots, we see that in the central part of the matter epoch (cf. at the same time the 
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Figure 2: All types of late-time solutions in the F = 1/B model as a phase diagram of the 
deceleration q and the dimensionless Hubble rate h = cH > 0, c = \2p\/ j3\^/'^ . The thick grey curve 
corresponding to q — Q goes through the unstable de Sitter fixed point aX q = —1 (circle) and the 
stable fixed point at q = ai ~ —0.74 (square). The arrows signify the direction of cosmic time. 

g-plot in the interval where q = 1/2) the EOS of the dark energy tends to stay around 
ojesiz) ~ 0, whereas deep in the radiation epoch {q = 1) it tends to values around ujes{z) ~ 
1/3. 

For the rest of this work, we want to avoid fractional powers of n in (4.5) and thus we 

4 _ 8 

take the value ?i = 3 implying Rcc Hi (x a s , which is still close to ACDM. 
5.3 The asymptotic late-time evolution 

In the previous sections we have analyzed the cosmic evolution in the matter and radiation 
epochs by making use of i? — )• 0, i.e. the approximate cancelation of both terms in (4.5). 
Accordingly, the behavior of pAcff was found without solving the complicated Einstein equa- 
tions (3.15). However, in the period from the recent (matter dominated) past to the asymp- 
totic future we cannot profit from the relation i? — )• 0, because the Hubble rate H is too 
small making the y-term in (4.5) completely negligible here, so that in good approximation 

B^'^-R^ + ]^g = 2AH\q - ^)(g - 2) . (5.9) 
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Obviously, the CC relaxation is supported now only by the smahness of H, since there are 
no terms in B which could cancel each other. Before we turn to numerical solutions of the 
Einstein equations in Sec. 5.4, let us discuss the very late-time evolution analytically, i.e. the 
asymptotic future regime {t — )• oo). For this purpose we consider first the model F = 1/B 
with y = - i.e. with B given by (5.9) - and neglect all energy density contributions except 
for pp and p\. This is a good approximation at late times since the Friedmann equation 
then reads 

3H^ 



Pm + PAcff = Pm + Pa + PF ^ Pa + PF 



(5.10) 



owing to Pm ~ — s- for t — t- oo. Moreover, the Hubble rate H must be very small 
in order to sustain the asymptotic late-time relaxation of the effective CC, and hence the 
effective CC tends to the asymptotically vanishing value of the critical density: pAcS = 
Pf + p\ ^ Pc ^ 0. In practice, this means that the equation pp + p\ = Q should be 
sufficient for analyzing the background evolution in the asymptotic regime. The induced 
term corresponding to the F = 1/B model can be computed explicitly with the help of 
(3.17). After a straightforward calculation, we find 

864 /3 H"^ 



PF 



IE", 



53 



5^2 - 3g - 5 



31q^ + 15q^ + 19q 



q (4g^ 



10) 

Wq + 7) 



q{Aq^ -Wq + 7))] 



(5.11) 



2i?4(2g2_5g_^2)2 2iJ5(2g2 _5g_^2)3 

For the sake of simplifying our discussion below, it is convenient to rewrite the previous 
equation as follows: 

ki q k2 

2H^ ~ Hb ' 2H^ 



PF = P 



(5.12) 



with 



ki 



5(g - ai){q - 02) 
ai,2 = ( 3 T V 



10 

%-6i)(g-62) 



109 



(5.13) 



bi,2 
2{q- 



Next, we express the asymptotic Friedmann equation p\ + Pf = in terms of a new 
dimensionless Hubble rate defined as h := cH and a scaled cosmic time r := t/c, where 



1 and 



respectively. Denoting q'{T) = dq/dr, we have q'{T) = cq{t) and 



hence Friedmann's equation can be cast as 

14 ki 



ko 



(5.14) 



62 -^'h^h' 

where "— " corresponds to x > and "+" to x < 0, respectively. Moreover, the time- 
evolution of h is the same as that of the true Hubble rate, 

h'{T) = -h^{q+l). (5.15) 
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From (5.13) and (5.14) it follows that g' and h' can be expressed as functions of only 
q and h, which fixes uniquely their evolutions independent of c. Instead of plotting a 
vector map, we show the integrated solutions in Fig. 2, where the arrows indicate the 
direction of time and the thick grey curve corresponds to q' = 0. In the following we will 
restrict our analysis to late time solutions (7 < ^ and H^h > 0, which implies 6, A;2 > 
and /ci ^ for g ^ oi —0.74, respectively. Consequently, the deceleration q decreases 
(becomes more negative) with time on the left-hand side and above of the (7' = curve and 
increases on the right-hand side and below it. Notice from (5.14) that the de Sitter regime 
q = -I satisfies li^ = ki/h'^ = 5(1 + ai)(l + a2)/81, and hence li = 27''^/^. As the Hubble 
rate decreases for q> —\ and increases for g < — 1, it follows that the de Sitter fixed point 
{q = —l,h = 27"^^^) - see the circle in Fig. 2 - is unstable, and a de Sitter final phase 
would require fine-tuning. On the other hand, the fixed point (g = ai ~ — 0.74,/i = 0) 
- see the square in the central part of Fig. 2 - corresponds to a stable final state similar 
to quintessence solutions. The red dashed curves describe universes with a decelerating 
past and a phantom (superaccelerated runaway) future, whereas the black curves approach 
the stable fixed point with quintessence-like behavior. Also the blue dashed-dotted curves 
end in the stable fixed point, however they are accelerating in the far past, possibly with a 
transient decelerating phase in between. Finally, the green dotted curves describe phantom 
universes with q < —1 all the time. 

Obviously the most interesting late time solutions obtained from the F{R,Q) = 1/B 
models under study are those showing stable asymptotic quintessence-like behavior. This 
kind of solutions also exist for more general models, as we shall see below, although the 
asymptotic future value of the deceleration q will be different in general. Let us look for 
accelerated solutions with constant g' > — 1 in the asymptotic future (i.e. with a terminal 
acceleration below the pure de Sitter regime). They correspond to power-law solutions of 
our field equations in the limit t — )■ 00, namely 

a{t) o^f{l + Ct'), (5.16) 

where a first-order correction \ Ct^\ <^ 1 (with ^ = const.) has been included^. We can search 
for these solutions within the generalized class of models of the form F(R, Q) = with 
m > 0. Obviously, the model we have been considering so far corresponds to the particular 
case m = 1. Direct calculations to first order in (^t^ via (3.26) and (3.17) lead to the 
following results: 

PF = ^ (PFO + PFI ■ Ct") , Pf = ^ (pfo + PFi ■ Ct") , (5.17) 
B = 12r^{9r^ - 9r + 2r)t-^ + 12rs(4 - 27r + 36r2 - 4s + Ors)*"^ • C^^ 

where {pfo-, Pfo) and (pfi, Pfi) are (dimensionless) time-independent terms corresponding 

more general approach for finding asymptotic solutions would be to start with the ansatz a{t) ~ 
t^{l + f{t)) and then linearize in the smooth function f{t). In this work, we will limit ourselves to explore 
the solutions of the form (5.16) as they are already able to capture the main features of the very late time 
behavior. 
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to the zeroth and first order corrections, respectively. In particular, 

K(r, m) = Or^ - 9r + 2 - 4:m^{9r - 4) + 3m{3r'^ - llr + 4). (5.19) 

However, these terms are in fact zero because in order to fulfill the late-time Einstein 
equation pp = —p\ we need K{r,m) = 0, otherwise pp ~ pFoB~^ oc t^™ would diverge 
for t — )• DO. In addition, we have to fix the correction term in the ansatz (5.16). This can be 
achieved by the choice s = —Am, as this insures that the term pp oc ppi • C,t'^B~'^ ~ ^s+4m 
remains essentially constant in time - and this also entails pp = const. Likewise, with this 
choice of s we have C,t^ = ^f"^*" — )■ with increasing t because m > 1, which means that 
the correction term in the solution (5.16) becomes smaller and smaller in the asymptotic 
regime, irrespective of the particular value of the coefficient C,. Accordingly, we find 

= r^'^'-Zl2f - - ^'^^ - + ('-'"^ 

-2m(-8 + 54r - llTr^ + Slr^)] , 

-2m2(-80 + 408r - 585r2 + 162r^) - 4m(-16 + 102r - 189r2 + 54r2 + 81r'^)] , 

implying 

16m2(4m - 1)(2 + 7n(4 - 9r) - 9r + 9r2) 
m+m = 3r2(9r2-9r + 2) ^^-''^ 

Since the zero-order terms ppQ and ppQ vanish, we obtain the correct asymptotic EOS for 
the F-terms, 

^P = PL = PEl = _i + 0(^^-4^) . (5.23) 

PF PFl 

It follows that the effective entity that we have called the "cosmon" , and which is respon- 
sible for the induced DE - of (modified) gravitational origin - behaves asymptotically as a 
cosmological term (up to very small corrections). Recall that this is actually so, although 
in a lesser extent, during most of the cosmological history prior our time. After all the 
cosmon "duty" is to continually neutralize the initial cosmological term leaving a small 
dynamical remainder - the measurable CC term. What we have just shown here is that, 
ultimately, it behaves (very approximately) as a true cosmological constant, with a value 
equal (but opposite in sign) to the initial CC term: pp{t — )• oo) = —p\- Finally, the 
condition ppQ oc K{r,m) = determines the leading power-law exponent r in the scale 
factor: 

r = (6m + 6)"' {I2m^ + 11m + 3 ± \/l44m4 + 200m3 + 81m2 + 10m + 1) . (5.24) 

Being r > 0, the evolution law (5.16) for the scale factor is an increasing one. For m = 1, 
we find the solutions r ~ 0.43 and r ~ 3.91, which are already close to the large m limit 
solutions r — 7- I and r — ?• 4m. Neglecting the correction term, we have a ~ r(r — 1)^"^, 
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and hence only the r > 1 solutions produce acceleration (a > 0), and only these can be 
interpreted as asymptotic quintessence-like solutions. 

Let us indeed consider in more detail the EOS of the effective vacuum energy /OAeff = 
Pf + Pa tti^ asymptotic regime. As we know, p^cs is essentially zero in this regime, 
which means that the corresponding EOS function cjgff cannot be easily derived from the 
original definition (3.34) because the latter involves the ratio between the two terms l-\-ujp 
and Pf + p\ both of which are very close to zero. Let us thus turn to equation (3.38) and 
apply it to the very late time epoch, where pm can be neglected. It is then easy to show 
that that equation can be recast in the form 

^^Ae. := ^ = , (5.25) 

Pc 1 - 3^^^) 

where Wg^'' is the asymptotic value of lOcs- At the same time, we know that in this regime 
PAcH = Pf + Pa tends to become arbitrarily close to pc (both going to zero with the same 
pace). Thus ^IacS — ^ 1, and consequently 

4P = -^ + lil + ^)- (5-26) 

This is the quintessence-like behavior of the EOS for the effective vacuum energy pAefr in 
the asymptotic limit. Notice that since the accelerated expansion implies coes < —1/3, we 
consistently obtain q < 0. Let us now compute explicitly the deceleration parameter q in 
the asymptotic regime. Using q = —1 — H /H^ and working out H = a/a directly from 
(5.16), we find 

H{t) = J (l - ^Ci"^"" + 0(a"^™)^) , (5.27) 
q{t) = -l + l + ^J<^_^(^t-^m ^ o{Qt-'^f . (5.28) 

We reconfirm from the last expression that q < because we look only for the solutions 
satisfying r > 1 - the time-dependent term in (5.28) also respects this feature because 
m > 1 and so it contributes negatively. Inserting the previous expression for q{t) in (5.26) 
we obtain the asymptotic EOS: 

J^){t) = -1 + f + ^^^^(l^ct-'" + OiCt-'-^f . (5.29) 
6r 6r'^ 

From this formula, it becomes apparent that the solutions r > 1 of equation (5.24) lead 
to quintessence-like behavior since u^s then lies in the interval — 1 < oj^s < —1/3 (with 
again a very small time dependence which does not alter this conclusion). Furthermore, 
if we take e.g. the case m = 1, equation (5.24) gives r ~ 3.91 for the r > 1 solution, and 
correspondingly equation (5.28) yields q ~ —0.74 and (5.29) provides uj^,^^ — —0.83 (up 
to very small time-dependent effects) - see figures 1 and 2. This value of UcS can be read 
directly from Fig. 1 - it corresponds to the value of tJeff (z) at the terminal point z = —1 
(i.e. t — oo). 
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We may finally derive the very late-time time behavior of pAeff- As we know, it behaves 
as pAes{t) Pcit) ^ ioi t ^ oo , but we can compute precisely how it decays with time 
in the last stages of its evolution. Substituting the asymptotic Hubble rate (5.27) into the 
asymptotic Friedmann's equation (5.10), we arrive at the desired result: 

PAeff = P\+PF= ( 1 - —Ct-'"' + 0{Ct-'n') ^ 7^ (t ^ OO) , 



SttGn \ r ^ ' J SttGn 

(5.30) 

where m > 1 for all the relaxation models under consideration. Let us also note (a posteri- 
ori) that the obtained result is consistent with our assumption that the asymptotic matter 
density decays always faster than the vacuum energy density, as we supposed from the very 
beginning in (5.10). Indeed, being p^ ~ cl~^ ~ t"^*^' - where we can neglect here the small 
correction term in (5.16) ~ the matter density is always much smaller than the effective 
vacuum energy pAcfT ~ because the minimum value of r (for accelerating solutions) is 
approximately r ~ 3.91, which is attained for m = 1 - cf . equation (5.24). For m > 1, r is 
larger than the previous value and the desired condition is always secured. 

Looking at the numerical examples presented in Figs. 1 and 5-8, we can check that 
the asymptotic t — )• oo (i.e. z — )• —1) EOS behavior uJcS = ^cs{z) of the effective CC is 
always quintessence-like (— 1 < Wcff < —1/3). However, in some cases (e.g. in Figs. 1 and 
5), it presents an effective phantom phase cjcfr < — 1 for non-asymptotic regions, actually 
for regions not far away in our past. Therefore, the accessible region to our observations 
could present an effective phantom-like DE behavior. In contradistinction to scalar models 
of the DE, this phantom-like behavior has nothing to do with negative kinetic terms, as 
it is of purely effective nature. For example, since Qcsiz) > in the affected regions of 
Figs. 1 and 5, the phantom behavior is just triggered by the fact that the set of terms 
in the field equations that we have collected under the name of the cosmon produces 
(1 + ojp{z))pf{z) < in those redshift segments ~ cf. equation (3.34). 

5.4 Parameters and initial conditions 

Our modified gravity action (3.7) with F = 1/B or more generally F = R^/B™' con- 
tains two parameters /? and y, where the latter is directly related to the radiation-matter 
transition. To solve the Einstein equations, we have to impose one more initial condition 
because F(R, Q) theories have more degrees of freedom than general relativity. In fact, 
while in the standard ACDM model the initial value of the deceleration is fixed simply 
by the current relative densities ri^, see equation (3.39) (specifically qo = ^2^/2 — Q\ for 
the present epoch), in the relaxation model we have to provide also the value of go- To see 
this, note that the induced F-term is a function of the type 

PF = PF(.H,q,q;/3,n,y) , (5.31) 

which depends in general on the three free parameters (4.7) and involves H, q and also q. 
Recall that n = 3 was fixed because it provides a sufficiently smooth transition from the 
radiation to the matter epochs, so only /3 and y remain as free parameters to be chosen. 
An explicit example is the late time form given in equation (5.11), where the q dependence 
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Figure 3: Possible future final states in the class of F = 1/B models. Shown is the deceleration 
parameter versus redshift for different values of the dimensionless parameter ^ = 3l3/{p\ Hq), 
assuming ^ < in all cases. Depending on the value of |(^|, in the asymptotic limit i — > oo (i.e. 
z — )> — 1) the universe expands approximately according to a power-law a ^ t*" (r ~ 3.91 for this class 
of models) or else it develops a future singularity characterized by a runaway acceleration (q <^ — !)■ 
Notice that de Sitter space-time {q — —1) lies on the boundary between the two possibilities. 



is manifest. Therefore, if we are given p\ and we fix the current value of the CC density, 
this amounts to fix the value of 

nO _ Pa + PF{Ho,qo,qo; ^,n,y) 

^^Aeff — -fi > {0.6Z) 

Pc 

which entails a non-trivial relation between the parameters. This relation involves explicitly 
the value of and in our case qo is fixed from equation (3.40). Therefore, equation (5.32) 
fixes qo, which is tantamount to say that qo must be provided as an additional input in 
order that the value of the current f^^gg is correctly matched. 

For definiteness, let us now focus for a while on the canonical type of models charac- 
terized hy F = 1/B, i.e. our starting class of models defined in section 3.2. As we know, 
for these models /3 has mass dimension 8 and can be represented as (3 = A4^, with A4 some 
mass scale. This suggests that the numerical analysis of the F = 1 / B models can be more 
conveniently performed using the dimensionless combination ^ = 6/(xHq) = 3/3/(/0^i?Q) 
as an input parameter. In Figs. 3 (^ < 0) and 4 (^ > 0), we plot the deceleration q for 
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Figure 4: As in Fig. 3, but for ^ > 0. Notice that here all curves have an asymptotic power-law 
expansion a(t) ^ f for t — > oo. However, not all curves are admissible since ^ is limited from below 
> 210) in order to have a matter era {q w i) in the past. 

different values and signs of ^ (or equivalently of qo). In these examples we use go — —0.6 
for the current deceleration, although one has some freedom here, too. The plots show 
a broad range of different solutions for q, with some of them close to the ACDM evolu- 
tion. Remember that the latter makes a precise prediction for the transition redshift from 
deceleration to acceleration, given by 



= -1 + ^2 ^- 0.75, (5.33) 

where the numerical value corresponds to our choice = 0.27 (hence fi^ = 0.73 for a flat 
universe). We see that most of the examples plotted in Fig. 3 - corresponding to different 
values of the parameter ^ and the initial conditions as in Fig. 1 - predict a transition 
redshift smaller and hence nearer to our time, meaning that the accelerated expansion is 
more recent. For large and negative however, one finds z* ~ 0.7 and the ACDM result 
can almost be matched - see e.g. the ^ = —10^ curve in Fig. 3. In this example, if we take 
as usual ~ with Mx ~ 10^'^ GeV, and we recall that Hq ~ lO'''^ GeV, we find 
M = ~ 10-4 eV, i.e. a mass scale characteristic of light neutrinos. Similarly in Fig. 4, 
corresponding to ^ > 0, although in this case we can see that there are values that give a 
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Figure 5: Model F R/B with n = 3, y = 0.7 x 10'^ , p\ = -lO^OGeV*, = 0.27, 

= 10^"*, qo w —0.6, qo — —0.5 Hq. The curves have the same meaning as in Fig. 1. 



transition redshift much earlier in time (i.e. a more remote onset of the acceleration period), 
although the scale Ai remains of the same order of magnitude. In fact, all values of |^| in the 
wide interval 1 ^ |C| ^ 10*^ render Ai in the light neutrino range 10~^ eV ^ A4 ^ 10~^ eV. 
However, for ^ > the values which are below 210 are clearly excluded as they correspond 
to (7 < in the past and therefore the matter epoch (q = 1/2) would not have occurred. 

A few more observations are worthwhile in regard to Figs. 3 and 4. For negative values 
of the asymptotic evolution follows either the power-law expansion discussed above 
(if I'^l ^ 105) a future singularity (if |^| < 104). This singularity corresponds to a 
situation leading to superaccelerated cosmology (q <^ —1) ending in a Big Rip. The 
boundary between both final states is de Sitter space-time {q = —1). Larger values of |.^| 
lead to an expansion behavior quite close to ACDM in the recent past. In the case of 
positive ^, only the power-law solution exists for t oo. However, as already mentioned, 
in this case |^| is bounded from below. Summarizing, late-time expansion histories close 
to ACDM (as suggested by observations) require large values of the parameter |^| (hence 
of /3, for a given p\), which always imply a power-law future state. In the limit |/3| — )• oo 
we find for both signs of p\ the value ~ — 0.57i/o) and the evolution of q becomes very 
insensitive with respect to /3, see Fig. 4. The last statement can be derived from (5.11), 
where we see that the terms in square brackets must vanish for |/3| — )• oo, otherwise the 
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Figure 6: Model F = R^/B with n = 3, y = 0.7 x IQ'^ "^'^ , p\ = -lO^^GeV'^, = 0.27, 
= 10"'', qo w —0.6, qo — —I.SHq. There is a future singularity at z > —1. The curves have the 
same meaning as in Fig. 1. 



Einstein equation (3.35) is not fulfilled. Consequently, q does not depend on (3 since it can 
be expressed function of q and H only. 

The phantom universes with q < —1 at all times (green dotted curves in Fig. 2) are not 
among the numerical solutions since we fixed the initial conditions to qo = —0.6. This also 
explains the absence of the solutions plotted as blue dashed-dotted curves in the case ^ < 0. 
Finally, let us also remark that although the numerical examples displayed in Figs. 1 and 
5-8 correspond to p\ < (basically because, as we have discussed in section 2, this is the 
situation in the important case of the SM of particle physics), a similar set of relaxation 
curves would be obtained for > 0. 



5.5 Generalized relaxation models: = W /B^ and beyond 

Without much effort the results in Sec. 4.2 can be carried over to many models which have 
a factor B^^ in pAcg and pAcflF- The natural generalization is to consider models of the 
form 

F':= = -p^ , (s > 0,m > 0; n > 2) . (5.34) 
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Figure 7: Model F ^ R^/B"^ with n = 3, y = 0.75 x lO^^ ^/a^ ^ _^Qeo QgV^ 17°, = 0.27, 
= 10"'', go ~ —0.6, qo — —O.SHq. The curves have the same meaning as in Fig. 1. 



Recall that we usually choose n = 3 to smooth the transition between radiation and 
matter eras. The canonical model with which we have started in section 3.2 is just the 
particular case F^, and the models considered in equation (5.17) are the F° ones. The 
cosmic evolution of q{z) in the matter and radiation epochs for does not show any 
relevant differences with respect to F^ because only the structure of B is important there, 
in the sense that the relaxation relation — )• holds. However, the relative energy 
densities ^n{z) of the matter /energy components and the EOS ujcs may evolve differently, 
which is demonstrated in Figs. 5, 6 and 7. Moreover, at late times also the deceleration q 
may differ significantly among the models. 

Apart from dynamical properties, the dimensionality and the value of the parameter /3 
depend on the particular generalized class F^ of F{R, Q) functional, too. One can easily 
show that the parameter /? is related to a power of a mass scale M as 

= M^-2^+^™. (5.35) 

All these models yield induced DE terms of the form pf,Pf with p > as a result of 

equations (3.17) and (3.18), where the derivatives of F and G introduce more factors of B 
in the denominators. One can readily show that the CC relaxation and tracking properties 
in the matter and radiation eras follow again from the relaxation condition. 
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Figure 8: The unified dark matter/dark energy scenario from Sec. 5.6 in the model F = /B^ 
with n = 3, J/ = 0.75 X IQ-^ H~^'^ , p\ = -1060GcV^ f)^ = 0.04, n° = lO"", sa -0.6, 
qo = —0.8Hq. The curves have the same meaning as in Fig. 1 except for fi^ (black dashed-dotted) , 
which represents the baryons only. The unified dark component exhibits a positive energy density 
PAcff > at all times. 

For estimating f3 let us assume that each time derivative in the expression for E^q - see 
equation (3.17) - yields a factor Hq today, and no accidental cancelations are present so that 
B ~ Hq. Consequently, in the model = l/B, we find ~ 13 /Hq for the current value 
of the induced term, as in the toy model example of section 4.1. Indeed, in the late 
epoch, where q ceases to play a significant role, the realistic model based on (4.5) produces 
qualitatively the same picture as in the simple one (4.2). Therefore, within one order of 
magnitude, equation (4.4) holds again and we can estimate the value of the parameter /3 
using the current value of the Hubble rate, Hq ~ 10~^^ eV. Since ~ —p\ and ~ 
we get 

|/3| =M^^ \p%\ ~ , (5.36) 

where /3 is expressed as a power of an energy scale Ai as in (3.12). Taking the standard 
GUT scale Mx ~ 10^^ GeV, we obtain once more a mass value for 7W of the order of a 
light neutrino mass: 

M ~ y^Mx Ho ~ 10^^ eV . (5.37) 
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Incidentally, this value for M is the geometric mean of the two most extreme mass scales 
available in our universe below the Planck mass. Such value is by the way not far from 
the mass scale associated to the current value of the CC: uia = {p%)^^^ ~ 10~^ eV, and it 
should therefore be considered quite natural. Actually, the value of Ai obtained in a given 
model could be larger than our rough estimation, especially if cancelations in E^q occur, 
but as we shall see there is a wide spectrum of possible values for M. in the class of 
models, and all of them within a reasonable particle physics range. 

In fact, applying the same arguments for the models = R/B and F2 = E? /B^ we 
find in both cases |/3| = ^ p\Hq Hq, which implies M ~ 0.1 GeV. This is again 

a reasonable mass scale in particle physics, as it is of the order of the characteristic QCD 
scale of the strong interactions, Ad ~ ^qcd — 100 MeV. Finally, the model F^ = B? / B 
yields |/3| = Ad^ = \pW ~ Af^' ^^^^ ^^^^ would be close to the initial vacuum 

energy scale Mx of the GUT ^ . This situation is also perfectly reasonable inasmuch as Mx 
is another natural scale of the problem, just the one placed at the starting point of the 
evolution. Remarkably, the mass parameter M of the relaxation mechanism is not only 
completely free from fine-tuning problems, it also lies in a perfectly reasonable range of 
particle physics masses, possibly related to neutrinos, QCD or even GUT models. In all 
these situations the relaxation mechanism takes care automatically that the observable CC 
scale of the vacuum energy density is pAeff(-ff) ^ p\ sA- H = Hq. 

In summary, the energy scale Ai of the CC relaxation mechanism could be very well 
related to reasonable mass scales in Particle Physics, thereby avoiding tiny energy scales 
Hq ~ 10^'^^ eV which so often appear in many dark energy models (for example, in 
quintessence models [13]). Once the correct F-term has been identified, the mass scale 
M should be considered a sort of "constant of Nature", and in this sense we expect it 
should bear direct relation with a particle physics scale of the SM or of a typical GUT 
scale. 

We remark that the fact that the initial cosmological term p\ was considered constant 
(but otherwise arbitrary), is not a real limitation for the efficiency of our dynamical ad- 
justment mechanism. This mechanism works for virtually any cosmological model in the 
general class of AXCDM models [24, 25]. In particular, we cannot exclude that the ini- 
tial vacuum energy density can be evolving with time through some cosmological quantity 
^ = '^(i), i-e. PA = pA{£,{t)). For example, it has been suggested in the literature that pA 
could be a quadratic function of the expansion rate: p\{H{t)) = uq + 112 H'^{t), where uq 
and 722 are both non-vanishing coefficients, see e.g. the recent papers [14, 15, 33] and the 
older ones [16, 17]. The reason for this "running" of the vacuum energy stems from the 
expanding background, and these papers suggested that the quantum effects provide pre- 
cisely this running behavior. For a summary of alternative proposals of evolving vacuum 
energy that have recently been tested in the light of the most recent cosmological data, see 
e.g. [22]. 

In the present framework, the effective vacuum energy pAcS is anyway a mildly evolving 

^Notice that although for this model the function / B itself does not increase with i/ — >■ (it remains 
constant in this limit), the derivatives of that function in equation (3.17) do indeed increase for — >■ and 
hence they are entirely responsible for the relaxation process in this case. 
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quantity, even in the absence of the aforementioned running contributions. But if one 
accepts that the running terms are also there, nothing essential is changed in our relaxation 
mechanism. One has to keep pA 7^ in the covariant conservation law (3.5) and the F- 
functional will adjust automatically the contribution of the additional terms. Technically, 
however, let us clarify that if these contributions are added directly at the level of the 
field equations, then in order to keep matter conservation we have to allow for a "running 
gravitational coupling Gn" as well. This possibility has also been contemplated in the 
literature and it corresponds to the so-called type-II AXCDM models [25] , in which the 
energy density conservation law involves pA 7^ and Gn ^ simultaneously. We refer the 
reader to the aforesaid reference for details. 

Finally there is still another modification of our relaxation functional that may be 
necessary in order to make our cosmological model a bit more realistic: we have to insure 
that the modified gravitational action (3.7), once it is written in a metric amenable for Solar 
System tests (namely, the Schwarzschild metric in the presence of a cosmological term), 
is able to pass these tests and at the same time keeps its ability to reduce dynamically 
the value of the cosmological constant within observations. While a detailed treatment of 
this problem will be presented elsewhere [57], in Appendix C we briefly describe the kind 
of modification that we have to introduce on the functional J-" such that it fulfills all these 
requirements without altering in any essential way the dynamical relaxation mechanism 
that we have studied so far. 

5.6 Special scenario: unified dark matter and dark energy 

In this section we discuss a special scenario of unified dark matter and dark energy. At least 
for the background evolution this seems to be possible because the dark energy density PAeff 
behaves like dust matter in the matter era, see equation (5.3). Therefore, it might replace 
dark matter completely. Whether the effective dark energy/matter in our setup provides 
the correct clustering properties to successfully seed structure formation will be studied in 
a future work together with its evolution on solar system and galactic scales^. If it works, 
we would find a MOND-like theory since our modified gravity model had to simulate the 
effect of standard dark matter. 

In the following, we concentrate on the cosmic background expansion. The model 
shown in Fig. 8 describes a universe with only baryons = 0.04, dark energy and 
radiation $7^ = 10~^. Since standard dark matter is absent, the relative energy densities of 
dark energy and the baryons sum up to unity in the matter era. At late times and during 
the radiation epoch there is no big difference between this special relaxation scenario and 
the previously studied relaxation models with dark matter. 

Finally, we mention two advantages of this approach. First, we do not need the pa- 
rameters related to dark matter (^^2^) because y fixes the radiation-matter transition. 
Second, it is possible to have pAcff > all the time, which leads to an effective EOS for 
PAcflf without spurious divergences. In fact, these divergences should not represent a real 
problem, they just indicate that E^q in (3.17) vanishes. However, in this special scenario 
the full 'g-component in the Einstein equation remains always positive. 

®In fact, there exists already some recent work considering the clustering effects of the vacuum energy [52] . 
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6. Conclusions 

In this work, we have made a thorough attempt to deal with the old cosmological constant 
(CC) problem [3] - the toughest cosmological conundrum of all times. It suggests that 
our vacuum should contain a huge energy density p\ coming from quantum zero-point 
energies and phase transitions in the early universe. This formidable problem cannot 
be solved by just finding a source for the late-time acceleration of our cosmos. Instead, 
it requires a powerful mechanism capable of disarming the large CC during all known 
stages of the cosmological evolution, not just at the last stages. Scalar field models, for 
instance, completely fail to solve the old CC problem, as they actually introduce two 
severe theoretical diseases: 1) extreme fine-tuning of the parameters (through e.g. an ad 
hoc counterterm in the potential that cancels the initial CC), and 2) an extremely tiny mass 
scale needed to match the ground state potential with the present value of the CC. Such 
mass scale is usually as small as -f^o ~ 10~^^ eV, i.e., very many orders of magnitude smaller 
than the natural mass scale m-A = (p^)^''^ ~ 10~^ eV associated to the observationally 
measured value of the vacuum energy density. Quintessence models, therefore, do not 
really solve the problem since they input two highly unnatural ingredients that completely 
spoil the credibility of the proposal. If that is not enough, they simply put the vacuum 
energy under the rug as if it would not exist, and then focus exclusively on the special 
properties of the newly invented scalar field. 

Clearly, a radically different new approach is required. In an attempt to make a first 
step in this direction, we have proposed a self-adapting (dynamical) relaxation mechanism 
based on a pure modification of gravity. Thanks to it, we can avoid fine-tuning at all 
stages of the cosmological evolution. More specifically, we have demonstrated that an 
arbitrarily large initial CC can be relaxed automatically by complementing the Einstein- 
Hilbert action with a class of action functionals of the Ricci scalar and the Gaufi-Bonnet 
invariant, J^{R,Q). The two most remarkable achievements of our relaxation mechanism 
is the complete absence of fine-tuning in the parameters of the model, and at the same 
time the absence of tiny mass scales. The modified gravity action induces a dynamical 
dark energy component, pp, associated to the J^-functional. This induced DE density acts 
effectively as a "dynamical counterterm", i.e. one that self-adjusts to the necessities of the 
universe at any given moment, and is able to efficiently neutralize any exceedingly large 
vacuum energy left over near the radiation epoch by the primeval inflation mechanism, 
thus preventing the universe from the disastrous effects that would ensue otherwise (as e.g. 
the disruption of the primordial nucleosynthesis epoch). Indeed, by becoming sufficiently 
large to compensate for the initial CC, the effective energy density pp generated by the 
J-"- functional self-adapts automatically to the initial conditions of the universe and to the 
subsequent radiation and matter epochs, and enforces the universe to follow the standard 
FLRW cosmic expansion history in each one of these epochs. 

As a result, the "residual" dark energy of our model (or effective vacuum energy 
density pAcs) appears, at every stage of the cosmological evolution, as a mildly evolving 
function of the expansion rate, PAcsiH) = p\ + pf{H). This function takes values much 
smaller than the initial p\ because there is a large (and automatic) compensation of it 
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with the values taken by pf{H), which have opposite sign. Moreover, |pAeff(-f^)| ^ \Pa\ 
holds good for all values of H at (and below) the radiation epoch, and for any given /9^. In 
addition, the current value of PAeff can be of the order of the measured one by the modern 
cosmological data, i.e. PAefT = PAesiH = Hq) ~ ~ 10~^''GeV^, without introducing 
fine-tuning or any unnaturally small parameter. In fact, the small PAefr the observable 
cosmological term in the present framework, whereas the large initial CC value (and the 
corresponding dynamical counterterm pp) cannot be resolved individually from the usual 
CC measurements. 

All in all, our model - the Relaxed Universe - appears very similar to the standard Big 
Bang model and resembles the late universe near our time (concordance ACDM model) 
in most respects, irrespective of the initial vacuum conditions. Therefore, the J-{R, Q)- 
cosmology appears all the time as a FLRW model with essentially constant and tiny cosmo- 
logical term. In addition, the vacuum energy pAeff exhibits remarkable tracking properties 
because it evolves like radiation or dust matter in the corresponding cosmological epochs. 
This feature of the model sheds considerable light also on the cosmic coincidence problem 
inasmuch as it provides a raison d'etre for the fact that the DE and DM densities are still 
close at present. Indeed, one finds that the DE and DM were cosmic companions with 
similar EOS and approximately proportional densities, and this was so all the time until 
very recently, namely until the matter epoch extinguished (meaning that q departed much 
from 1/2) and the universe - which is constantly fueled with the big initial p\ - was forced 
to use the last resort of the relaxation mechanism available under these circumstances, 
which is to take a very low value of H in order to compensate for p^ and still keep show- 
ing a small observable PAeff- That crucial event decoupled once and forever the tracking 
behaviors of matter /radiation and DE, and since then the DE dominates the universe ex- 
pansion. Such breakdown of the tracking property at the end of the matter epoch, on the 
other hand, pushed the subdominant DE density curve upwards until crossing the decaying 
DM one at some point after the matter epoch and hence near our time. This is the origin 
of the "coincidence" event in our late time neighborhood. But, most remarkable of all, 
the primary trigger of this crossing event (representing the startup of the DE era) is the 
relaxation mechanism itself, which is therefore simultaneously responsible for explaining - 
and linking together - the two fundamental CC conundrums of our current Cosmos: the 
old CC problem and the cosmic coincidence problem. 

Worth noticing is also the fact that the asymptotic late time behavior of the universe 
in the DE epoch mimics the quintessential or phantom behavior. Only under special initial 
conditions, the future behavior can be de Sitter space-time, but any perturbation of these 
conditions would tilt its evolution into one of the aforementioned forms. 

The dynamical solution found in this paper resembles, in a sense, the one so longed for 
by scalar field model builders, except that these kind of models always failed because they 
ultimately involve fine-tuning (as proven by the famous Weinberg's no-go theorem [3] ) . 
Even if the J-{R^ t/)-cosmology lacks at the moment of a fundamental understanding of 
the origin of the J^-functional, the situation is not very different from the innumerable 
and unsuccessful attempts based on scalar field models where one has no fundamental 
motivation for using a particular potential that serves our purposes. The advantage of 
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the J-"(-R, ^)-cosmology, in contrast, is that the job is done by gravity itself, and is done 
correctly, namely without fine-tuning inconsistencies. This was never achieved with scalar 
field models nor with any other cosmological model that we know of. 

In its current stage, we have been able to devise a relaxation mechanism which starts 
working automatically in the early radiation dominated epoch, hence right after the reheat- 
ing process that follows primordial inflation. It was not our purpose here to identify the 
fundamental physics that may have triggered this mechanism. It could have been caused by 
an "effective inflaton" reaching the minimum of its potential, especially if the minimum lies 
below zero and having an arbitrary value |yo\|. This value fixes the initial (usually large) 
CC of the universe just before starting the radiation epoch. The latter case appears much 
more likely than just assuming (as usually done in the literature) that the inflaton potential 
has been fine-tuned to almost vanishing ground state energy at the radiation epoch. Fur- 
thermore, since our relaxation J-"- functional is constructed from a modified gravity theory, 
it is natural to complement it with the standard B? terms of the renormalizable effective 
action of gravity. These terms do not alter the relaxation mechanism in any significant 
way, as we have seen, but they can just take care of the UV behavior of the theory and 
even be responsible for inflation itself. 

In the framework of modified gravity, the successful relaxation of the CC requires an 
action functional -F(i?, Q) with a special structure, namely one which leads to an induced 
dark energy density being able to compensate dynamically the initial CC. However, our 
approach is not limited to a very specific model of this kind; actually, a whole class of 
working relaxation models has been uncovered. On the basis of several analytical and 
exact numerical examples, we have shown that they differ mostly at late times in accordance 
with the requirement of the current cosmos being close to ACDM. Interestingly enough, our 
action functional J- contains a mass scale M whose natural value can either be a typical 
particle physics mass of the standard model of strong and electroweak interactions (say, a 
neutrino mass, the QCD scale etc) or even a typical Grand Unified mass Mx- But M is 
never an extremely tiny mass scale of order of Hq ~ 10"'^'^ eV, which so often appears in 
the DE models of the literature. In summary, we can say that the typical values for the 
fundamental mass scale M that are needed to solve the CC fine tuning problem in our 
framework range from the present CC mass scale, m\ = (p^)^''^ ~ 10"'^ eV, to the value 
of the initial CC mass scale, {p\f'^ Mx^ lO^^ GeV. What else could be more natural? 

The model universe we have envisaged here falls within the general class of the AXCDM 
models of the cosmic evolution [24, 25, 26], in which a new entity X (the "cosmon") appears 
in interplay with the cosmological term. Such an entity is in general not a field, but a 
complicated effective quantity. The cosmon is in fact associated here to the induced DE 
density pp which is generated from the J-"- functional at the level of the field equations. 
While in Ref. [24] the cosmon served the main purpose of solving (or highly palliating) 
the cosmic coincidence problem, here the cosmon is eventually responsible for solving the 
fine-tuning CC problem, and hence it provides a first fundamental step towards solving 
(or highly alleviating) the old CC problem. Amazingly enough, it provides a clue for 
solving the coincidence problem as well, as we have seen. Therefore, while in the present 
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realization the cosmon does not appear as an elementary particle, it does however fully 
accomplish its main aim, which is to cure the CC fine-tuning problem. Let us recall that 
the name "cosmon" was coined for the first time in the literature in Ref. [51], and in its 
first implementation it was a scalar field, hence afflicted by the aforementioned Weinberg's 
no-go theorem. In the modern version, it is not. Therefore, the new cosmon well deserves 
its name, for it fully adapts to the original sense of the old literature as being an entity 
which is able to dynamically adjust the value of the cosmological constant [51]. The cosmon, 
indeed, generates the dynamical counterterm pp which is able to continuously neutralize 
the very large ~ in all epochs of the post-inflationary cosmic evolution. 

In summary, in this paper we have extensively discussed the working principle of the 
relaxation mechanism and the corresponding background evolution of the J-{R^ Q) cosmol- 
ogy. In doing this we have found an exceptional model universe, the Relaxed Universe^ in 
which all the main tensions with the observational data are relaxed and can be smoothly 
accommodated in its inner dynamics. The Relaxed Universe seems to be a promising 
candidate for our real universe since it is essentially free from the toughest cosmological 
puzzles, the old CC problem and the cosmic coincidence problem. Its explicit construction 
constitutes a "proof of existence" of a dynamical mechanism that may account for the value 
of the current cosmological constant, starting from a value of arbitrary size. However, al- 
though our results are very encouraging, this is just one more step in the right direction. 
More efforts are necessary to uncover all the main properties of our model and resolve 
upcoming problems. In the future, we plan to investigate further aspects of this setup, 
e.g. the evolution of perturbations, the Newtonian limit and possible instabilities, leading 
to new insights and improvements. Also, it will be interesting to know if the unification 
of dark matter and dark energy can be fully realized in this framework. And of course 
one should eventually look for a fundamental theory of the underlying action functional, 
which at the moment is lacking. 

The previous issues are certainly difficult problems of fundamental nature, which must 
be addressed in the future in order to assess whether the model (or some variation of it) 
can eventually become a realistic one. In the meanwhile, the very theoretical construction 
and analysis of the Relaxed Universe does show that it is possible to envision a model 
universe in which an arbitrarily large vacuum energy early deposited in it by quantum 
field theory or string theory dynamics can be rendered innocuous in a fully dynamical 
way - hence without any fine-tuning - with the pure work of gravity, and still resulting in 
an universe whose cosmological features (in particular the tiny value of the total vacuum 
energy) resembles to a great extent the standard ACDM model. In this sense, and in spite 
of the many difficulties that may lie ahead us, one could say that this model universe should 
be a first step towards solving the old CC problem. 
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A. F{R,Q) modified gravity 

In this section we derive the equations of motion of the f{R, S, T) and F{R, Q) modified 
gravity theories in the metric formahsm. In the modified action the / and F functionals 
contain the Ricci scalar R = g°-^Rab as well as S" = RabR""^ and T = RatcdR"'^'^'^ , which 
are scalar invariants built from the Ricci Rab and Riemann Rabcd tensors. Q is the Gaufi- 
Bonnet invariant Q = R? — AS + T. In general, f{R, S, T) theories will introduce new 
degrees of freedom, which potentially lead to instabilities not present in general relativity. 
However, some problems can be avoided by specialising to F[R, Q) functionals involving 
only the Ricci scalar and the Gaufi-Bonnet invariant. For instance, all functionals of the 
form f{R,T — 45) are ghost free [47], and in particular this is obviously the case of the 
F[R,Q) functionals. With respect to the more general f{R,S,T) ansatz this is not a 
strong restriction on an FLRW background, since in terms of the Hubble rate H and the 
deceleration q, the invariants have the form 

R = 6H^{l-q), S = 12H^{q^ -q+1), T = 12H^ {q^ + I) , ^ = -24/f^g. 

Therefore, S and T can be replaced by = — and T* = — Q, respectively. 
Note that these replacements may not hold on the level of the equations of motion or in 
general metrics. 

First, we fix the notation. Our metric Qmn has the signature (+,—,—,—), and the 
Riemann and respectively the Ricci tensors are given by 

-r>a -pa -pa _i_ pa pm pa pm t-» jya 

^ bed bd,e be,d ' mc-*- bd md^ bci ^bc -ft bca> 

where F^^ are the Christoffel symbols. In short, our conventions here are characterized 
by the three basic signs (— ,+, — ) relative to the sign conventions of Ref. [53] concerning 
metric, Riemann tensor and the r.h.s. of Einstein's equations respectively, following the 
well-known classification scheme proposed in that reference. 

For a general f{R, S, T) modified gravity theory the action reads 

S = j (fx^\f{R,S,T), 

and its variation 5S with respect to the metric yields 

5S[f{R,S,T)] = j d^x [5i^\)fiR,S,T) + ^\iF'''6R + F'^dS + F^5T) 

= J d^xy^\[-^gabf{R,S,T) + RabF'^ + F%-gabaF'^ 

+ 2F^RaeRb' - 9ab{F^R'"');ed + 2{F^ R'%ae9bd " a{F^ Rab) 

+ 2F^RarstRr' - 4{F^Rm;mn]Sg''' + (surface terms), (A.l) 
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where S{y^\) = -^y^lgabSg""^ and 

From SS[f{R, S, T)] one can easily derive the variation 6S[F{R, Q)] after applying in (A.l) 
the replacements 

f{R, S, T) ^ F{R, G), F^^ + F^ = 2RF^ + F^, 

^ = -4F^, F^ ^ I^F^ = F^, (A.2) 

where on the righthand side of the arrows F^'^ correspond to the partial derivatives of 
F(R,Q). On a FLRW background the results of this procedure are given in Sec. 3. 

B. The CC fine-tuning problem in QFT: a more detailed account 

In quantum field theory the CC fine-tuning problem (discussed only at the classical level in 
section 2) becomes much harder. We feel that despite the widespread criticisms scattered 
over the literature against the fine-tuning procedure, chiefly in connection with the CC 
problem, a detailed discussion of it is lacking, or at least is not usually made available to 
the general reader (apart from some vague notions deprived of a minimum quantitative 
analysis), and therefore here it seems the right place to pay due attention to it^. Let us 
flesh out the new ingredients of the problem at the quantum level in this Appendix. As a 
first important step, we need to renormalize the theory because otherwise it is impossible 
to reach finite results. While it is unavoidable to have a curved background in the presence 
of a vacuum energy, we shall confine our renormalization discussion to flat space in order 
to avoid a cumbersome presentation. Even without curvature, it will become apparent 
that the presence of quantum effects makes the fine-tuning problem extremely difficult 
and demands a cleverer solution to the problem. We start by recalling some basic issues 
concerning the renormalization of the effective potential [54] . Indeed, in QFT the field f 
becomes a quantum field operator ip and therefore its ground state value < f > must now 
be interpreted as the vacuum expectation value (VEV) of this operator: < ip >= (0|(^|0). 
Notwithstanding, the theory can still be handled as if were a classical field provided its 
potential (2.1) is renormalized into the effective potential, V — >• V^s, where 

v^^ = v + hVi + h'^V2 + h^V3 + ... (B.i) 

The quantum effects to all orders of perturbation theory arrange themselves in the form 
of a loopwise expansion where the number of loops is tracked by the powers of h. Thus, 



^Let us clarify that the fine-tuning problem is not at all privative of the CC approach to the dark 
energy. In most papers dealing with the DE from the point of view of quintessence, the fine-tuning problem 
is usually ignored ab initio as if this problem would not go with them. This is however an unfortunate 
misconception. These models are plagued, too, with fine-tuning problems and in fact in no lesser extent 
than the traditional CC approach - including other niceties, such as the unwanted presence of extremely 
tiny mass scales. 
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at one loop we have only one power of h, at two loops we have two powers of h etc. For 
h = 0, however, there are no loops and the effective potential just reduces to the classical 
potential, V, given by equation (2.1) in the electroweak standard model. On the other 
hand, the loop terms in (B.l) can be split into two independent contributions, one having 
loops with no external legs (vacuum-to- vacuum parts Vp^) and the other having loops 
with external legs of the quantum matter field if (i.e. the loop corrections V^^l^i^) to the 
classical potential): 

V, = Vi^^ + V^2M , V, = vP + kS(v.) , V, = + V^ii^)... . (B.2) 

As a result, the effective potential (B.l) at the quantum level splits naturally into two 
parts, one which is independent and another that is (/3-dependent: 

V,s{^) = ^zPE + V,,^i{ip) , (B.3) 

where 

^zpE = h vi^^ + v^^ + h^vP + .... (B .4) 

is the zero-point energy (ZPE) contribution, which consists in the sum of all the vacuum- 
to- vacuum parts of the effective potential. The ZPE part is sourced exclusively from closed 
loops of matter fields (i.e. vacuum loops without external (/?-legs). For instance, at one-loop 
order, only the first term on the r.h.s. of (B.4) contributes. In dimensional regularization, 
it gives 

r(l) _ * , 4-n / " . r 7,2 I 21 _ ^ ,,4-n ' " 



^ 2^ J (27r)" L J 2'^ 7 (27r)"-i ^ 

where we have performed a Wick rotation {dko = idk^) into Euclidean space. Notice that ^ 
is the characteristic 't Hooft mass unit of dimensional regularization, and we used n — )• 4 in 
the final result. The second equality in (B.5) - in which we performed the contour integral 
on /c4 - was only to make more transparent the connection of this expression with the more 
traditional form of the ZPE of the vacuum fluctuations. Finally, the expression 



4 



is the one-loop /3-function of the vacuum term [16]. As warned from the very beginning, 
these results are valid only in flat space-time, and are given just to illustrate in a concrete 
way the structure of the expansion (B.4). In the presence of a curved background, the 
expression (B.5) must be supplemented with contributions from the curvature invariants. 
The effect of the matter fields on them is to make their coefficients to run as effective 
charges^. In all the cases the vacuum-to- vacuum part does not depend on Lp, but only on 



*For a discussion of these curvature terms and tlie potential implications on the vacuum energy in an 
expanding background, see [14] . 
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the set of parameters P = m, A, ... of the classical potential. The ZPE receives in general 
contributions to all orders of perturbation theory, except at zero loop level since Vzpe is a 
pure quantum effect that vanishes for h = 0. Besides, there is the (/^-dependent part: 



This one is not purely quantum (i.e. it does not vanish for ?i = 0) as the first term is not 
proportional to ^ - it corresponds to the tree-level contribution V{ip) or classical potential. 
In the electroweak standard model, it is given by (2.1). The above 99-dependent part of 
Ves receives in general also contributions to all orders of perturbation theory, and vanishes 
for = since in this case all the loops have external (p legs, including the tree-level 
part. Thus, Vesi^^p = 0) = Vzpe, which is a number. In other words, (B.7) constitutes 
the quantum corrected effective potential (excluding the ZPE number). However, the full 
effective potential contains both contributions. 

In the above discussion, all the field theoretical ingredients (m. A, (p and Ves) are in 
reality bare quantities {niQ, Aq, (fo and T4fro) that require renormalization inasmuch as the 
loopwise expansion is UV-divergent order by order. However, renormalization just means 
that we replace all the bare quantities with renormalized ones (in some given renormal- 
ization scheme with a specific set of renormalization conditions) plus counterterms (which 
are also scheme dependent and are partially fixed by the condition of cancelling the UV- 

1 /2 

divergences): mo = m + 6m, Aq = A + 6X, (po = ip = {1 + 6Z^/2) ip... Of course, a 
similar splitting occurs with the vacuum term, which was originally a bare term pAvacO- 
This parameter is of special significance in the present discussion. We must also split it 
into a renormalized piece plus a counterterm: pAvacO = PA vac + ^PAvac- The full set of 
counterterms is essential to enable the loop expansion to be finite order by order in per- 
turbation theory. For instance, if we would renormalize the theory in the MS scheme in 
dimensional regularization, the suitable counterterm reads: 



The next crucial point is to apply the basic renormalization recipe in QFT, which says 
that we must equate the full bare effective action (EA) to the full renormalized EA after 
performing the aforementioned renormalization transformation. Since the effective poten- 
tial is obtained from the effective action in the limit of constant mean field or background 
scalar field [38, 54], it follows that the bare EA boils down to 



where O is the total space-time volume. After performing the renormalization transforma- 
tion of parameters and fields, we obtain r[Lp{fi),m{fi), X{fi), p\{fj,), fi]. This renormalized 
EA must be the same as the bare one. In particular, in the constant mean field limit 
the two terms on the r.h.s. of (B.9) are replaced by pAvac(M) and Vcs[^p{fi),m{fi), X{n), n] 



i^scai(v^) = v{^) + hv};>,{^) + v;S(v.) + v}^{^) + ... 



(B.7) 






(B.9) 
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respectively, where dpAvac has been additively incorporated in the structure of the latter. 
Since the two overall expressions must be equal, we have 



Mvaco + ycs{vo,mo, Ao) = PAvac(Ai) + Kfr(v^(/w), rn{ii), A(//); fi) . (B.IO) 

Notice that the renormalized result depends on an arbitrary mass scale fi (say, the formerly 
used 't Hooft mass unit in dimensional regularization, or an arbitrary momentum p in a 
momentum subtraction scheme). Furthermore, in the above expressions, the counterterms 
involve some regulator (e.g. space-time dimensionality n in dimensional regularization, or a 
simple cutoff) and can be chosen to cancel all the divergences produced by the computation 
of the loop corrections to the potential. The regulator, therefore, eventually disappears, 
but the mass scale fi remains as a reflex of the arbitrariness of the subtraction point. For 
instance, in the free theory or in the absence of SSB the one-loop MS-renormalized vacuum 
energy density in flat space is obtained from (B.4) and (B.IO) for ip = 0: 

PA = PAvacO + VzPEO = PAvs,c{p) + VzPE(pt) , (B-H) 

where 

Vzpe{p) = hvj:^^ +5 pA,,c (B.12) 

is the one-loop renormalized ZPE, which, as we see, necessarily requires the counterterm 
associated to the vacuum term to become a finite quantity. From (B.5) and (B.8) the 
renormalized vacuum energy density for f = can be easily derived: 

PA = PAvM + ^ln-^--j . (B.13) 

We see that it is explicitly dependent on p, but this dependence represents only an internal 
parameterization of its structure, which is helpful to track the various kinds of quantum 
effects [14]; the overall /U-dependence, however, must eventually cancel. In fact, the renor- 
malized parameters are finite quantities which are also functions of p: ip = (p{p), m = m{p), 
A = X{p), /OA vac = PAva.c{p), and since the vacuum energy cannot depend on the arbitrary 
scale p, the sum of the renormalized vacuum term and the renormalized potential must 
be globally scale-independent (i.e. //-independent). This is obviously so because the bare 
vacuum term and bare effective potential were scale-independent to start with. Thus, from 
(B.IO) we have 

p-j-[PAviLcim{p),X{p);p)) + Vcs{p{p);m{p),X{p);p)] =0. (B.14) 
dp 

This relation implies that the full effective potential is actually not renormalization group 
(RG) invariant (contrary to some inaccurate statements in the literature), but it becomes 
so only after we add up to it the renormalized CC vacuum part pAvac- In reality, the 
structure of the effective potential (B.3) is such that the previous relation splits into two 
independent RG equations: 

P^[PA^;a.c{m{p), \{p); p)) + VzpE{rn{p), X{p); p)] =0 (B.15) 
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and 

/"^ Kical('/3(/i);m(/u), A(;u);^) = 0. (B.16) 

Equation (B.15) shows that it is only the strict vacuum-to-vacuum part (i.e. the ZPE) 
the one that needs the renormahzed vacuum term pAvac to form a finite and RG-invariant 
expression, whereas the renormahzed (/9-dependent part of the potential (i.e. the tree-level 
plus the loop expansion with external y^-tails) is finite and RG-invariant by itself. This is 
of course the essential message from the renormalization group ~ see Ref. [14, 19] for an 
appropriate physical interpretation. Explicitly, equation (B.16) reads 

d d d ^ 



1 12 

where as usual Pp = fidP/d^ [P = m, A, ...) and = ^dlnZ^ /(9/i. Similarly, equation 
(B.15) can be put in the form (B.17), except that the ip term is absent. 

The following interesting result now emerges. Plugging equation (B.13) in the general 
RG equation (B.15), we find immediately that the renormahzed vacuum term /OAvac(/^) 
must obey the following one-loop RG-equation: 

rfpAvac _ m^h _ (1) 

which is of course the reason why we called the expression f5^^ in (B.6) the one-loop /3- 
function of the vacuum term. Notice two particular things of the one-loop result: A is not 
involved, and m does not run with /x. Both things cease to be true at higher orders. 

The one- loop renormalization of the effective potential is standard [54] , although the 
usual discussions on this subject rarely pay much attention to disentangle the ZPE part 
from it. Let us do it. Once more we equate the bare and renormahzed effective poten- 
tials after having performed the renormalization transformation of parameters and fields. 
Sticking to the MS scheme in dimensional regularization to fix the counterterms, one finds: 

where from (2.1), 

V"{^) = m^ + ^X(p'^ . (B.20) 
This expression already contains the renormahzed ZPE, as it is just its value at (/? = 0: 

Notice that this result is perfectly consistent with (B.13). As remarked, the full effective po- 
tential is not RG-invariant, in particular its one-loop form (B.19). The truly RG-invariant 
expression is the sum /?Avac(/^) + V^fr(v')) as indicated by (B.14). The vacuum energy den- 
sity is just the VEV of this expression. Therefore, we are now ready for addressing the 
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CC fine-tuning problem in the context of a well-defined, renormalized, and RG-invariant 
vacuum energy density in flat space. 

Once the full effective potential has been renormalized, the two loopwise expansions 
(B.4) and (B.7) become finite. Furthermore, the basic equation (2.14) remains formally 
the same in the quantum theory, i.e. the physical energy density associated to the CC is 
the sum of the vacuum part plus the induced part. The only difference is that the induced 
part now contains all the quantum effects, i.e. it reads pAind = (K;ffo('/5)), or equivalently 
PAind = {^cS^if))^ where V^g^{ip) = 14ff ((/?(//); m(//), A(/x); //) is the renormalized effective 
potential. Notice that the latter includes the (renormalized) ZPE part, which was absent 
in the classical theory. Thus, the physical CC emerging from the renormalization program 
(in any given subtraction scheme) reads 

PAph = PAvac + PAind = P^c + (^^(V')) = PA^ac + V^^E + (Kc'al(^)) , (B.22) 

where for simplicity we have obviated the p-dependence in the renormalized parts. The 
first equality in (B.22) can be considered the bare result and is formally identical to the 
last equality after removing the superscript "ren" everywhere. Although we have provided 
explicit one-loop results for the sake of concreteness, the above formulae are completely 
general. They also illustrate the systematics of the renormalized perturbative expansion 
of the vacuum energy and the kind of different contributions that are generated to all 
orders of perturbation theory, not just to one-loop order. In this sense, they are sufficiently 
descriptive to formulate now the severity of the CC problem and its relation with fine- 
tuning. 

Indeed, since the expression (B.22) is the precise QFT prediction of the physical value 
of the vacuum energy to all orders of perturbation theory, it follows that it must be equal 
to the observational measure value [1, 2], i.e. pAph = Pa - 2-8 x IQ-^^GeV^ We have 
already seen in section 2 that the lowest order contribution from the Higgs potential is 55 
orders of magnitude larger than p^, and that this enforces us to choose the vacuum term 
PAvac with a precision of 55 decimal places such that the sum pAvac + PAind gives a number 
of order lO^^^GeV^. The problem is that the fine-tuning game, ugly enough as it is at 
the classical level, becomes actually much more perverse at the quantum level. Indeed, 
recall that we have the all order expansions (B.4) and (B.7). Therefore, the quantity that 
must be equated to p^ is not just (2.12) but the full r.h.s. of (B.22), which is a finite 
and RG-invariant expression. In other words, instead of just requiring to satisfy equation 
(2.15), we must fulfill the much more severe one: 

10-^' GeV^ = PAvac - 10^ GeV^ + ^il^'^ + h"" vj,^^ + vj,^\... 

Clearly, since on the r.h.s. of this equation there are two independent perturbatively 
renormalized series contributing to the observed value of the vacuum energy (viz. the ZPE 
series and the series associated to the quantum corrected Higgs potential), the selected 
numerical value for the "renormalized counterterm" pAvac must be changed order by order 
in perturbation theory; specifically, pAvac must be re-tuned with 55 digits of precision as 
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many times as indicated by the order of the highest loop diagram providing a contribution 
to the CC that is smaller than the experimental number on the l.h.s. of equation (B.23). 
Let us assume that each electroweak loop contributes on average a factor (/^/(IGvr^) times 
the fourth power of the electroweak scale v = (if) ~ 100 GeV (see section 2), where g 
is either the SU (2) gauge coupling constant or the Higgs self-coupling, or a combination 
of both. It follows that the order, n, of the highest loop diagram that may contribute 
to the measured value of the vacuum energy, and that therefore could still be subject to 
fine-tuning, can be approximately derived from 



Take e.g. g equal to just the SU{2) gauge coupling constant of the electroweak SM, which 
satisfies g"^ / (IGvr^) = Oem/ (^tt sin^ Ow) — 2.5x 10^"^. This is actually a conservative assump- 
tion because in practice there are larger contributions in the SM associated to the big top 
quark Yukawa coupling, but it will suffice to illustrate the situation. Since v = (y?) ~ 100 
GeV, we find n ~ 21. Therefore the 21th electroweak higher order loop diagrams still 
might contribute sizeably to the value of the CC, and must therefore be readjusted by an 
appropriate choice of the renormalized value of the vacuum term /OAvac- This is of course 
preposterous and completely unacceptable. Even though this is nonsense, it is implicitly 
accepted by everyone that admits that such "technical trick" is a viable solution to the CC 
problem. It goes without saying that this situation worsens even more for higher energy 
extensions of the SM, such as in grand unified theories (GUT's). 

We see that the presence of the bare vacuum term /OAvac was not optional. While we 
could have dispensed with it in the classical theory, it is no longer possible in the quantum 
theory. This term is essential in order to provide both a UV-divergent counterterm to 
renormalize the ZPE part of the effective potential, as well as a finite (and fine-tuned) 
counterterm to compensate the large induced CC contribution. In other words, the vacuum 
term is not just a classical object but an indispensable ingredient insuring the technical 
consistency of the theory at the quantum level. However, even if technically possible, this 
fine-tuning game does not seem to be a terribly convincing recipe to solve the CC problem 
in any sensible way. Specially in the QFT case, where we have made evident that it becomes 
a nightmare of incommensurable proportions; in fact one which cannot be sustained on any 
sound theoretically ground. 

We hope to have convinced the reader that a dynamical solution of the old CC problem 
is absolutely mandatory, especially if we do take the vacuum energy of QFT any seriously. 
Of course, one can hide the vacuum energy of the SM of the strong and electroweak inter- 
actions under the rug and resort to a particular scalar field and claim that the cosmological 
constant adopts its value entirely from its ground state. However, this way of proceeding 
does not explain, or give any hint of why we can "happily" ignore the huge vacuum energy 
density predicted by the most successful quantum field theory at our disposal, viz. the Stan- 
dard Model of the Elementary Particles, nor does explain at all why the value of the energy 
density of such field is precisely of the order of 10"^'' GeV^ rather than, say, Mj^ ~ 10^^ 
GeV^, which is what everyone would expect for a non-SM field connected with Grand Uni- 




(B.24) 
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fied Theories. Therefore, from our point of view, it is much more convincing to envisage an 
exphcit mechanism deahng with the entire value of the vacuum energy (irrespective of its 
size) and reducing it (dynamically) to its present small value, even if the fundamental origin 
of this mechanism is unknown. After all, in spite of the tremendous variety of DE models 
(including all the scalar field models of course) , there is not a single one which can be shown 
to emerge from a fundamental theory, say GUT's or strings, and none of them provides 
a solution to the fine-tuning problem either. Despite scalar fields (dilaton-like ones, for 
example) are usually associated with string theory, they have a priori no serious impact 
on the CC problem without providing (by hand) some suitable form of the corresponding 
scalar field potential, and subsequently exerting (also by hand) a strong fine-tuning of its 
parameters. In contrast, the Relaxed Universe presented in this paper provides at least an 
efficient dynamical scape to fine-tuning, irrespective of the size of the initial vacuum energy 
in the early post-inflationary universe. Thus, hopefully, this prototype model should be a 
first indispensable step towards eventually finding a fundamental theory of the cosmological 
constant. 



C. Solar System environment 

So far we have shown that the CC relaxation mechanism works well in a cosmological 
setup. In the following we briefly discuss how the mechanism operates in the Solar System 
environment, which we describe by the Schwarzschild-de Sitter (SdS) metric in spherical 
coordinates, i.e the standard Schwarzschild metric including the presence of a cosmological 
term A: 

ds'^ = Adt^ - A-^ dr^ - r^dO^ - W 9 dcj? , (C.l) 
A = \-'±-\kr^^\-'^-^, (C.2) 

where the Schwarzschild radius of a body of mass m is denoted by = 2Gm, and r^s is 
the characteristic size of this finite universe (for a given value of A). Indeed, for = 
(or in practice for r ^ r^) the above line element corresponds to a 3-dimensional sphere 
of radius r^g = y^S/A. Obviously, for the observed value of A ~ Hq this is a very large 
radius of the order of the Hubble horizon, Vfis « Hq^, and for r much smaller than this 
value the presence of the term r'^/r'^g in the metric can be neglected. Therefore, we have 
to insure that in our relaxation model the effective A remains within its observable value 
A ~ Hq also in the above SdS metric. 

On the background metric (C.l) the Ricci scalar, the Gaufi-Bonnet term and respec- 
tively the denominator function B = ^R"^ + are given by 

R=12r-^, g = 2Ar-l + I2rlr-^, B = 108 r'^^ + Gr^r'^ (C.3) 

One can easily see that those terms in (C.3) that do not depend on the presence of the mass 
m just follow from the corresponding cosmological results (3.20) and (5.9) upon setting the 
correspondences H — ?■ r~^g and q — t- —1, as expected from the behavior in the de Sitter 
space limit. 
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In order to obtain the metric (C.l) as a solution of the generahzed Einstein's equations 
(3.15) of our framework, we minimally extend our action functional by adding a new term 
responsible only for the CC relaxation in the Solar System. For instance, consider the 
model 

^-M' (C.4) 

where the j3/B part has been discussed before with /? ~ Pa^q given in Eq. (5.36). With 
the help of Eq. (A.l) it is possible to determine the extra contribution to the tensor Eat 
in (3.15) following from the term /3q/R in J-" computed in the metric (C.l). One finds 
^ab ~ 9ab{~ Po'^'ds / •> ^"^^ ^'^^ show bclow that f3/B induces a negligible correction in 
the Solar System environment. In the Einstein equations (3.15) we safely neglect all terms 
except and p\ in accordance to our procedure in Sec. 5.3. This implies 2E'^^+gahp\ = 0, 
with all corrections suppressed by the large CC term p\. Consequently, we find I5q ~ 
p^r'^g ~ P\Hq for the Solar System parameter, which allows the comparison of both terms 
in T. As a result, the condition \fi/B\ <^ \I3q/R\ holds for radii r smaller than 

the critical radius = (^^5^*)^^^' which is sufficiently large due to r^s ~ . Thus, 
neglecting j3/B m. T \s justified in this environment, and the CC is successfully relaxed 
by the new term only. In contrast, in the cosmological background the newly added term 
I3q/R ^ P\{Hq/H)'^ is much smaller than the first term of (C.4) in the matter era because 
we have /3/S ~ /j^ as a result of the relaxation mechanism discussed in Sec. 4. Furthermore, 
let us also remark that the late time behavior will be controlled once more by the original 
term /3/i? ~ in the functional, due to its higher power of H in the denominator. 

It means that the cosmological evolution of the modified model (C.4) will be essentially 
identical to the original one for the recent past, present and future. 

Finally, let us discuss the potential occurrence of fifth forces coming from the scalar 
degree of freedom in the J-' = Pq/R action functional. Since these forces are part of 
the gravitational sector, they have to show up in the metric that follows from solving the 
Einstein equations. In our case, we have seen above that the Schwarzschild-de Sitter metric 
in (C.l) is a good approximate solution. Fifth forces can only emerge from the corrections 
to that metric, which we found to be strongly suppressed in the Solar System environment. 
Note, that our gravitational action with a large CC term p\ differs significantly from 
the model J- (x {R — p^/R) from Refs. [55, 56], where the late-time acceleration follows 
from an interplay of both terms in T in the absence of the CC. In the latter setup, the 
parameter p. ~ Hq is very small and the mass of the scalar degree of freedom is of the same 
order of magnitude. Hence, it will mediate a long-range fifth force which is in contrast to 
observations [58]. In our scenario, the parameter corresponding to p \s /3q, which is much 
larger than in the aforementioned situation (in fact, (3q oc p\Hq) owing to the relaxation 
condition operating on the big initial /j^. One can show that this also leads to a much 
larger value for the physical scalar mass. For p\ in the ballpark of a typical GUT, the range 
of the interaction would be very short [57]. Of course, these arguments are consistent with 
Eq. (C.l) as a good solution in the Solar System. 

In summary, the simple extension of the modified gravity action by (Sq/R does not 
alter our previous results at the cosmological level, but it represents a possible way to relax 
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the CC also in the Solar System environment without introducing long-range fifth forces. 
A more detailed exposition of the last point will be presented elsewhere [57]. 
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